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ABSTRACT 

This paper focuses on cosmological constraints derived from analysis of WMAP data alone. 
A simple ACDM cosmological model fits the five-year WMAP temperature and polarization 
data. The basic parameters of the model are consistent with the three- year data and now better 
constrained: fif,/!^ ^ 0.02273 ± 0.00062, n^h'^ = 0.1099 ± 0.0062, = 0.742 ± 0.030, = 
0.963j:°:gJ^, T = 0.087 ± 0.017, and erg = 0.796 ± 0.036, with h = 0.719t!^:Jj2?- With five years of 
polarization data, we have measured the optical depth to reionization, r > 0, at 5ct significance. 
The redshift of an instantaneous reionization is constrained to be Zrcion — 11-0 ± 1.4 with 68% 
confidence. The 2a lower limit is Zicion > 8.2, and the 3(T limit is Zrcion > 6.7. This excludes 
a sudden reionization of the universe at z = 6 at more than 3.5ct significance, suggesting that 
reionization was an extended process. Using two methods for polarized foreground cleaning we 
get consistent estimates for the optical depth, indicating an error due to foreground treatment of 
T ~ 0.01. This cosmological model also fits small-scale CMB data, and a range of astronomical 
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data measuring the expansion rate and clustering of matter in the universe. We find evidence 
for the first time in the CMB power spectrum for a non-zero cosmic neutrino background, or 
a background of relativistic species, with the standard three light neutrino species preferred 
over the best-fit ACDM model with N^s = at > 99.5% confidence, and iVoff > 2.3 (95% CL) 
when varied. The five-year WMAP data improve the upper limit on the tensor-to-scalar ratio, 
r < 0.43 (95% CL), for power-law models, and halve the limit on r for models with a running 
index, r < 0.58 (95% CL). With longer integration we find no evidence for a running spectral 
index, with dus/dhik — — 0.037 ± 0.028, and find improved limits on isocurvature fluctuations. 
The current WMAP-only limit on the sum of the neutrino masses is ^ rrii, < 1.3 eV (95% CL), 
which is robust, to within 10%, to a varying tensor amplitude, running spectral index or dark 
energy equation of state. 

Subject headings: cosmic microwave background, cosmology: observations, polarization, early 
universe 



Introduction 



The Wilkinson Microwave Anisotropy Probe (WMAP), launched in 2001, has mapped out the Cosmic 
Microwave Background with unprecedented accuracy over the whole sky. Its observations have led to the 
establishment of a simple concordance cosmological model for the contents and evolution of the universe, 
consistent with virtually all other astronomical measurements. The WMAP first-year and three-year data 
have allowed us to place strong constraints on the parameters describing the ACDM model, a flat universe 
filled with baryons, cold dark matter, neutrinos, and a cosmological constant, with initial fluctuations de- 
scribed b};_nBarlY_^cale4m^M law fluctuations, as well as placing limits on extensions t o this simple 



model (ISpergel et al 



20031 . l2007f ). With all-sky measurements of the polarization anisotropy (jKogut et al 



20031 : iPage et al.ll2007r ). two orders of magnitude smaller than the intensity fluctuations, WMAP has not only 
given us an additional picture of the universe as it transitioned from ionized to neutral at redshift z ~ 1100, 
but also an observation of the later reionization of the universe by the first stars. 

In this paper we present cosmological constraints from WMAP alone, for both the ACDM model and 
a set of possible extensions. We also consider the consistency of WMA P constraints w i th oth er recent 
astronomical observations. This is on e of seven f ive-ye ar WMAP papers. iHinshaw et al.l (|2008l ) describe 
the data p r ocessi ng and basic results, iHill et al.l (j2008l ) present new be am models an d window functions. 
Gold et al.l ( 20081 ) describe the emission from Galactic foregrounds, and I Wright et al. the emission 



Komatsu et al.l (|200i 
data. 



from extra-Galactic p oint sources. The angular power spectra are described in Nolta et al. ( 2008h . and 



present and interpret cosmological constraints based on combining WMAP with other 



WMAP observations ar e used to produce ful l-sky maps of the CMB in five frequency bands centered at 
23, 33, 41, 61, and 94 GHz ( Hinshaw et al. 20081 ). With five years of data, we are now able to place better 
limits on the ACDM model, as well as to move beyond it to test the composition of the universe, details 
of reionization, sub-dominant components, characteristics of inflation, and primordial fluctuations. We have 
more than doubled the amount o f polarized data us ed for cosmological analysis, allowing a better measure 



of the large-scale E-mode signal ( Nolta et al. 2008h . To this end we test two alternative ways to remove 



Galactic foregrounds from low resolution polarization maps, marginalizing over Galactic emission, providing 
a cross-check of our results. With longer integration we also better probe the second and third acoustic peaks 
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in the temperature angular powe r spectrum, and have many more year-to-year difference maps available for 
cross-checking systematic effects (jHinshaw et al.ll2008l) . 



The paper is structured as follows. In 31] we focus on the CMB likelihood and parameter estimation 
methodology. We describe two method used to clean the polarization maps , describe a fa s t met hod for 
computing the large-scale temperature likelihood, based on work described in lWandelt et al.l (|2004[ ). which 
also uses Gibbs sampling, and outline more efficient techniques for sampling cosmological parameters. In fJ3] 
we present cosmological parameter results from five years of WMAP data for the ACDM model, and discuss 
their consistency with recent astronomical observations. Finally we consider constraints from WMAP alone 
on a set of extended cosmological models in 21 s-i^d conclude in ^ 



2. Likelihood and parameter estimation methodology 
2.1. LikeUhood 



The WMAP likelihood function takes the same format as for the three-year release, and software im- 
plementation is available on LAMBDA (http://lambda.gsfc.nasa.gov) as a standalone package. It takes in 
theoretical CMB temperature (TT), E-mode polarization (EE), B-mode polarization (BB), and temperature- 
polarization cross-correlation (TE) power spectra for a given cosmological model. It returns the sum of 
various likelihood components: low-£ temperature, low-£ TE/EE/BB polarization, high-£ temperature, high- 
£ TE cross-correlation, and additional terms due to uncertainty in the WMAP beam determination, and 
possible error in the extra-galactic point source removal. There is also now an additional option to compute 
the TB and EB likelihood. We describe the method used for computing the low-^ polarization likelihood 
in Section I2.1.1[ based on maps cleaned by two different methods. The main improvement in the five-year 
analysis is the implementation of a faster Gibbs sampling method for computing the £ < 32 TT likelihood, 
which we describe in Section r2. 1.21 

For ^ > 32, the TT likelihood uses the combined pseudo-Ci spectrum and covariance matrix described 



Hinshaw et al.l (|2007l ). estimated using V and W ban ds. We do not use the EE or BB power spectra at 



£ > 23, but continue to use the TE likelihood described in lPage et al.l (|2007l ) , estimated using Q and V bands. 
The errors due to beam and p oint sources are treated the same as in the three-year analysis, de s cribed in 
Appendix A of Hinshaw et al.l (|2007i ). A discussion of this treatment can be found in Nolta et al. ( 20081 ). 



2.1.1. Low-£ polarization likelihood 



We continue to evaluat e the exact likelihood for the polarization maps at low multipole, £ < 23, as 
described in Appendix D of iPage et al.l (j2007l ). The input maps and inverse co variance ma t rix us ed in the 
main analysis are produced by co- adding the template-cleaned maps described in Gold et al. ( 2008h. In both 
cases these are weighted to account for the P06 mask using the method described in lPage et al.l (j2007l ) . In the 
three-year analysis we conservatively used only the Q and V ban ds in the likelihood. W e are now confident 
that Ka band is cleaned sufficiently for inclusion in analyses (see iHinshaw et al.l (|2008l ) for a discussion) . 



We also cross-check the polarization likel ihood by using polari zation maps obtained with an alternative 



component separation method, described in Dunkley et al. ( 2008h . The low resolution polarization maps 
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in the K, Ka, Q, and V bands, degraded to HEALPix A^sidc = 8 0, are used to estimate a single set of 
marginalized Q and U CMB maps and associated inverse covariance matrix, that can then be used as inputs 
for the < 23 likelihood. This is done by estimating the joint posterior distribution for the amplitudes and 
spectral indices of the synchrotron, dust, and CMB Q and U components, using a Markov Chain Monte 
Carlo method. One then marginalizes over the foreground amplitudes and spectral indices to estimate the 
CMB component. A benefit of this method is that errors due to both instrument noise and foreground 
uncertainty are accounted for in the marginalized CMB covariance matrix. 



2.1.2. Low-£ temperature likelihood 



For a given set of cosmological parameters with theoretical power spectrum Cg, the likelihood function 
returns p{d\Ce), the likelihood of the observed map d, or its transformed aim coefficients. Or iginally, the 



likeli hood code was written as a hybrid combination of a normal and lognormal distribu tion (I Verde et al 
2003). This algori t hm did not properly r nodel the tails of the likelihood at low multipoles ( Efstathiou 20041 : 



Slosar et all l2004t lO'Dwver et al.l 120041: iHinshaw et all 120071) . and so for the three-year data the ^ < 30 
likelihood was computed exactly, using 



p(d|Q) 



exp[(-(l/2)d'^C-id] 
VdetC 



(1) 



where C is the c ovariance matri x of the data i nclu ding bot h the signal cova riance matrix and noise C(C^) = 
S(Q)+N fe.g.. lTegmari (|l99i : lBond et al.l (fligS) : Hinshaw et all (|2007l) V This approach is computation- 
ally intensive however, since it requires the inversion of a large covariance matrix each time the likelihood is 
called. 



In IJewell et all (|2Q04| ): IWandelt et all (|2004t ): lEriksen et all (|2Q04al ) a faster method was developed 
and implemented to compute p{d\Ce), which we now adopt. It is described in detail in those papers, so 
we only briefly outline the method here. The method uses Gibbs sampling to first sample from the joint 
posterior distribution p(Cf,s|d), where Ci is the power spectrum and s is the true sky signal. From these 
samples, a Blackwell-Rao (BR) estimator provides a continuous approximation to p{Ce\d). When a flat prior, 
p{Ci) = const, is used in the sampling, we have p{Ci\d) oc p{d\Ci), where the constant of proportionality 
is indep endent of C/. The BR estimator can then be used as an accurate representation of the likelihood. 



p{d\Cg) (jWandelt et al 



2004 



Chu et al.ll2005h . 



The first step requires drawing samples from p(Ce, s\d). We cannot draw samples from the joint dis- 
tribution directly, but we can from the two conditional distributions p(s|C^,d) and p{Ci\s,d), each a slice 
through the {Np x ^max)-dimensional space. Samples are drawn alternately, forming a Markov Chain of 
points by the Gibbs algorithm. For the case of one Ci parameter and one s parameter the sampling goes as 
follows. We start from some arbitrary point (C^,s*) in the parameter space, and then draw 



{C, 



i+l J+1' 



(2) 



by first drawing C^"*"^ fromp(Cf |s*, d) and then drawing s*+^ from p(s|C^^^, d). Then we iterate many times. 
The result is a Markov chain whose stationary distribution is p{Ci, s\d). This is extended to vectors for Cg 



^The number of pixels is 12N'^^^^, where A^side = 2'' for resolution 3 l lGorski et aLlbOOsh . 
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(of length ^niax) and s (of length Np) following the same method: 



p(s»|Q,d) 
p(Q|s'+\d), 



(3) 
(4) 



sampling the j o int di s tribution ip(C^ s|d). This samphng procedure is also described in I Jewell et alj (|2004l ): 
Wandelt et all |20o4 ): lEriksen et all |2007h . 



The first conditional distribution is a multivariate Gaussian with mean S*(S* + N)^^d and variance 
[(S*)^^ +N^^]^^, so at each step a n ew signal vector s'"*"^, o f size Nj,, can b e draw n. This is computationally 
demanding however, as described in lEriksen et alJ (|2004ar ); I Wandelt et alJ (|2004l ). requiring the solution of 
a linear system of equations Mv — w, with M = 1 + S^/^N~^S^/^. These are solved at each step using the 
conjugate gradient technique, which is sped up by finding an approximate inverse for M, a preconditioner. 
This requires computation of the inverse noise matrix, in spherical harmonic space, which is done by 
computing the components of term by term using spheric al harmonics in pixel space. There are more 
efficient ways to compute N^^ (jHivon et al.ll2002t lEriksen et a l. 2004b), but computing the preconditioner 
is a one-time expense, and it is only necessary to compute harmonics up to £ = 30. 

The second conditional distribution, p{Ce\s, d) is an inver se Gamma distr i bution , from which a new Ci 
vector of size £max can be rapidly drawn using the method in lWandelt et al.l (|2004l ). Sampling from these 
two conditional distributions is continued in turn until convergence, at which point the sample accurately 
represents the underlying distribution. This is checked in practice using a jacknife test that compares 
likelihoods from two different samples. 

Finally, once the joint distribution p{Ce, s|d) has been pre-computed, the likelihood for any given model 
Ci is obtained by marginalizing over the signal s, p{d\Ce) oc / p{Ce,s\d)ds, which holds for a uniform prior 
distribution p{C'e). In practice this is computed using the Blackwell-Rao estimator. 



p{d\Ce) cx / piCi\s)p{s\d)ds « — Vp(Q|s') 



(5) 



where the sum is over all uq samples in the Gibbs chain. Since p{Ce\s'^) ~ p{Cf\(j\), where — {2t 



\Sin 



and Sim are the spherical harmonic coefficients of s, one only needs to store al at each 



step in the Gibbs sampling. Then, each time the likelihood is called for a new Ci, one computes C = 
Yl^=iP{C'i\o'\) /no- This requires only 0(^rnax?^G) computations, compared to the full 0{Np) evaluation of 
equation [TJ This speed-up also means that the exact likelihood can be used to higher resolution than is 
feasible with the full evaluation, providing a more accurate estimation. 



Code details: Choice of £ limits, smoothing, and resolution 



The co de used for WMAP is adapted from the MAGIC Gibbs code described i nlWandel t' ('2003'):I Wandelt et al, 
(|2004l ). The input temperature map is the five-year ILC map described in iGold et al.i (|2008il . To produce 
correct results, the Gibbs sampler requires an accurate model of the data. This means that the signal covari- 
ance matrix S(C^) cannot be approximated to be zero except for multipoles where the smoothing makes the 
signal much less than the noise. For the full WMAP data set, this would require sampling out to ^ ~ 1000, 
with i Vsirifl = 512. This is co mputationally expensive, taking more than of order 10'^ processor-hours to con- 
verge (jO'Dwver et al.ll2004l ). Instead we reduce the resolution and smooth the data to substantially reduce 
the required multipole range, speeding up the computation. 
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Fig. 1. — This figure compares the low-^ TT power spectrum computed with two different techniques. At 
each I value, we plot the maximum likelihood value (tic mark), the region where the likelihood is greater 
than 50% of the peak value (thick line) and the region where the likelihood is greater than 95% of the peak 
value (thin line). The black lines (left side of each pair) are estimated by Gibbs sampling using the ILC 
map smoothed with a 5 degree Gaussian beam (at HEALPIX A^sidc — 32). The light blue line (right side of 
the pair) is estimated with a pixel-based likelihood code with A'sido = 16. The slight differences between the 
points are primarily due to differences in resolution. At each multipole, the likelihood is sampled by fixing 
the other Ci values at a fiducial spectrum (red) . 



The ILC map is smoothed to 5.0 degree FWHM, and sampled at N^idc = 32. The process of smoothing 
the data has the side effect of correlating the noise. Properly modeling the correlated noise slows down the 
Gibbs sampling, as it takes longer to draw a sample from p(s|Cf,d). We therefore add uncorrelated white 
noise to the map such that it dominates over the smoothed noise. However, the added noise must not be 
so lar ge that it changes th e likelihood of the low-£ modes; cosmic variance must remain dominant over the 



noise (jEriksen et al.ll2007n . so we add 2 /iK of noise per pixel. In Appendix A the noise power spectra are 
shown. 

The Gibbs sampler converges much more slowly in regions of low signal to noise. Because of this, we 
only sample spectra out to £ = 51 and fix the spectrum for 51 < £ < 96 to a fiducial value, and set it to zero 
for I > 96. The BR estimator is only used up to ^ = 32 for cosmological analysis, so we marginalize over the 
32 < £ < 51 spectrum. The likelihood is therefore p{L\d) = J p{L, H\d)dH, where L and H refer to the low, 
£ < 32, and higher, 32 < £ < 51, parts of the power spectrum. Examination of the BR estimator shows it to 
have a smooth distribution. Tests of the results to various input choices, including choice of resolution, are 
given in Appendix A. We note that by using the low-^ likelihood only up to £ < 32, this breaks the likelihood 
into a low and high £ part, which introduces a small error by ignoring the correlation between these two 
parts of the spectrum. However, this error is small, and it is unfeasible, in a realistic sampling time, to use 
the BR estimator over the entire £ range probed by WMAP. 
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Parameter Description 



ui), Baryon density, 

uic Cold dark matter density, Qch^ 

Dark energy density, with w = —I unless stated 
Amplitude of curvature perturbations at fcp = 0.002/Mpc 

Us Scalar spectral index at ko = 0.002/Mpc 

T Reionization optical depth 

Agz SZ marginalization factor 

drts/dlnk Running in scalar spectral index 

r Ratio of the amplitude of tensor fluctuations to scalar fluctuations 

a_l Fraction of anti-correlated CDM isocurvature (See Sec l4.1.3t 

ao Fraction of uncorrelated CDM isocurvature (See Sec l4.1.3)l 

Ngff Effective number of relativistic species (assumed neutrinos) 

uji. Massive neutrino density, Q^h^ 

f2fc Spatial curvature, 1 — Qtot 

w Dark energy equation of state, w = pu e / pD E 

Yp Primordial Helium fraction 

Xe Ionization fraction of first step in two-step reionization 

Zr Reionization redshift of first step in two-step reionization 

as Linear theory amplitude of matter fluctuations on 8 Mpc scales 

Ho Hubble expansion factor (lOOh Mpc^^km s^^) 

^rrii, Total neutrino mass (eV) ^m,j = 9iflu ti^ 

Qm Matter energy density fli, + Qc + 

fimh^ Matter energy density 

to Age of the universe (billions of years) 

Zreion Redshift of instantaneous reionization 

»710 Ratio of baryon to photon number densities, 10^^ {ni, / n-y) = 273.9 



Table 1: Cosmological parameters used in the analysis, http://lambda.gsfc.nasa.gov lists the marginalized values for 



these parameters for all of the models discussed in this paper. 

In Figure [T] we show slices through the Cg distribution obtained from the BR estimator, compared to 
the pixel likelihood code. The estimated spectra agree well. Some small discrepancies are due to the pixel 
code using iVgidc = 16, compared to the higher resolution A^sido — 32 used for the Gibbs code. 



2.2. Parameter Estimation 



We use the Markov Chain Monte Carlo methodology described in lSpergel et al.l (|2003l . 120071 ) to explore 
the probability distributions for various cosmological models, using the five-year WMAP data and other 
cosmological data sets. We map out the full distribution for each cosmological model, for a given data 
set or data combination, and quote parameter results using the means and 68% confidence limits of the 
marginalized distributions, with 



M 

d^xL{d\x)p{x)x, = -Y,^l (6) 
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where is the jth value of the ith parameter in the chain. We also give 95% upper or lower limits when 
the distribution is one-tailed. We have made a number of changes in the five-year analysis, outlined here 
and in Appendix B. 



We parameterize our basic ACDM cosmological model in terms of the following parameters: 

{uJb, t^c, f^A, T, Us, A^}, {Asz} 



(7) 



defined in Table [T] is the amplitude of curvature perturbations, and the spectral tilt, both defined at 
pivot scale kg = 0.002/Mpc. In this simplest model we assume 'instantaneous' reionization of the universe, 
with optical depth r, in which the universe transitions from being neutral to fully ionized during a change 
in redshift of z = 0.5. The contents of the Universe, assuming a flat geometry, are quantified by the 
baryon density Wf,, the cold dark matter (CDM) density ujc and a cosmo logical constant, f^A- We tre at the 
contribution t o the power spectrum by Sunyaev-Zeldovich fluctua tions (Sunvaev fc Zel'dovich 197Cll ) as in 
Spergel et al.l (|2007l ). adding the predicted template spectrum from lKomatsu fc Seliatg (|2002[ ). multiplied by 
an amplitude Asz, to the total spectrum. This template spectrum is scaled with frequency according to 
the known SZ frequency dependence. We limit < Asz < 2 and impose unbounded uniform priors on the 
remaining six parameters. 



We also consider extensions to this model, parameterized by 

{dus/dln fc, r, a_i, ofQ, il_R-, w, u^, N^s, >p, Xe, Zr) 



(8) 



also defined in Table [T] These include cosmologies in which the primordial perturbations have a running 
scalar spectral index dns/dlnk, a tensor contribution with tensor-to-scalar ratio r, or an anti-correlated 
or uncorrelated isocurvature component, quantified by a-i, uq. They also include models with a curved 
geometry flk, a constant dark energy equation of state w, and those with massive neutrinos '^rrii, = 
9An^h^ eV, varying numbers of relativistic species Ncb, and varying primordial Hehum fra ction Yp. There 
are also models with non-instantaneous 'two-step' reionization as in ISpergel et al.l (j2007[ ). with an initial 
ionized step at Zr with ionized fraction Xe, followed by a second step at z = 7 to a fully ionized universe. 

These parameters all take uniform priors, and are all sampled directly, but we bound A'cff < 10, w > 
—2.5, Zr < 30 and impose positivity priors on r, a_i, ao, Wi/, Yp, and as well as requiring < a;e < 1. 
The tensor spectral index is fixed at fit = — r/8. We place a prior on the Hubble constant of 20 < Hq < 100, 
but this only affects non-flat models. Other parameters, including as, the redshift of reionization, Zroion, 
and the age of the universe, to, are derived from these primary parameters and described in Table [T] A 
more extensive set of derived parameters are provided on the LAMBDA website. In this paper we assume 
that the primordial fluctuations are Gaussian, and do not consider constraints on parameters describing 
deviations from Gaus sianity in the data. These are discussed in the cosmological interpretation presented in 
iKomatsu et al-lboosl ). 

W e continue to use the CAMB code jlewis et al.lboooh . based on the CMBFAST code |Seliak fc Zaldarriaga 
1996 ). to generate the CMB power spectra for a given set of cosmological parameters H. Given the improve- 
ment in the WMAP data, we have determined that distortions to t he spec t ra du e to weak gravita tional lensing 
should now be included (see e.g. Isianchard fc SchneideJ |l987h : ISeljaM (|l997h : Eu fc Okamoto (2002). We 
use the lensing option in CAMB which roughly doubles the time taken to generate a model, compared to 
the unlensed case. 



^We use the pre-March 2008 version of CAMB that treats reionization as 'instantaneous' (with width z ~ 0.5) and fully 
ionizes hydrogen but not helium at reionization. 
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Fig. 2. — The temperature angular power spectrum corresponding to the WMAP-oniy best-fit ACDM model. 
The grey dots are the unbinned data; the black data points are binned data with la error bars including 
both noise and cosmic variance computed for the best-fit model. 



We have made a number of changes in the p arameter sampling m ethodology. Our main pipeline now 
uses an MCMC code originally developed for use in lBucher et al.l (|2004l ). which has been adapte d for WMAP. 



For in creased speed and reliability, it incorporates two changes in the methodology described in lSpergel et al 



(120071). It uses a modified sampling method that generates a single chain for each model (instead of the four, 
or eight, commonly used in cosmological analy ses). We also use an alternative spectral convergence test 
that ca n be run on a s ingle chain, developed in lDunklev et al.l (|2005l ). instead of the Gelman & Rubin test 



used in ISpergel et al.l (120071). These are both d escribed in Appendix B. We also use the publicly available 
CosmoMC sampling code ( Lewis &: Bridlel2002 ) as a secondary pipeline, used as an independent cross-check 
for a limited set of models. 



3. The ACDM Cosmological Model 

3.1. WMAP five-year parameters 

The ACDM model, described by just six parameters, is stil l an excellent fit t o the WMAP data. The 

With more observation 



temperature and polarization angular power spectra are shown in lNolta et al 



(200. 



the errors on the third acoustic peak in the temperature angular power spectrum have been reduced. The 
TE cross-correlation spectrum has also improved, with a better measurement of the second anti-correlation 
at £ 500. The low-£ si gnal in the EE spectrum, due to reionizatio n of the universe ( Ng &: Ng IQQSt 

The best-fit 6 



Zaldarriaga fc SeliakI 119971 ). is now measured with higher significance (jNolta et al.ll2008f ). 

parameter model, shown in Figure is successful in fitting three TT acoustic peaks, three TE cross- 
correlation maxim a/minima, and the low-^ EE signal. The model is compared to the polarization data in 
Nolta et al. ( 2008[) . The consistency of both the temperature and polarization signals with ACDM continues 
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Parameter 


3 Year Mean 


5 Year Mean 


5 Year Max Like 




2.229 ±0.073 


2.273 ±0.062 


2.27 




0.1054 ±0.0078 


0.1099 ±0.0062 


0.108 




0.759 ±0.034 


0.742 ±0.030 


0.751 


ris 


0.958 ±0.016 




0.961 


T 


0.089 ±0.030 


0.087 ±0.017 


0.089 




(2.35 ±0.13) X 10-3 


(2.41 ±0.11) X 10-3 


2.41 xlO-3 


as 


0.761 ±0.049 


0.796 ±0.036 


0.787 




0.241 ±0.034 


0.258 ±0.030 


0.249 




0.128 ±0.008 


0.1326 ±0.0063 


0.131 


Ho 


73.2+3;i 


71.9t|? 


72.4 


^rcion 


11.0± 2.6 


11.0± 1.4 


11.2 


to 


13.73 ±0.16 


13.69 ±0.13 


13.7 



Table 2: ACDM model parameters and 68% confidence intervals from the five-year WMAP data alone. The three- 
year values are shown for comparison. For best estimates of parameters, the marginalized 'Mean' values should be 
used. The 'Max Like' values correspond to the single model giving the highest likelihood. 



to validate the model. 



The five-year marginaliz ed distributions for ACDM, shown in Table[2]and Figures [3] and |4l are consistent 
with the three-year results (jSpergel et al.ll2007n . but the uncertainties are all reduced, significantly so for 
certain parameters. With longer integration of the large-scale polarization anisotropy, there has been a 
significant improvement in the measurement of the optical depth to reionization. There is now a 5(t detection 
of r, with mean value r = 0.087±0.017. This can be compared to the three-year measure of t = 0.089±0.03. 
The central value is little altered with two more years of integration, and the inclusion of the Ka band data, 
but the limits have almost halved. This measurement, and its implications, are discussed in Sec 13.1.11 

The higher acoustic peaks in the TT and TE power spectra also provide more information about 
the ACDM model. Longer integration has resulted in a better measure of the height and position of the 
third peak. The highest multipoles have a slightly higher mean value relative to the first peak, compared 
to the three-year data. This can be attributed partly to improved beam modeling, and partly to longer 
integration time reducing the noise. The third peak position constrains l^^^^^fe, while th e third peak height 
strongly constrains the matter density, flmh^ ( Hu fc Whitel 19961 Hu Sz Sugiyama 1995). In this region of 
the spectrum, the WMAP data are noise-dominated so that the errors on the angular power spectrum shrink 
as 1/t. The uncertainty on the matter density has dropped from 12% in the first year data to 8% in the three 
year data and now 6% in the five year data. The CDM density constraints are compared to three-year limits 
in Figure [3l The spectral index still has a mean value 2.5cr less than unity, with = 0.963^g'oj^5. This 
continues to indicate the pref erence of a red spectrum consistent with the simplest inflationary scenarios 
(^indc 200i iBovle et al.ll2006f ). and our confidence will be enhanced with more integration time. 



Both the large scale EE spectrum and the small scale TT spectrum contribute to an improved measure 
of the amplitude of matter fiuctuations. With the CMB we measure the amplitude of curvature fluctuations, 
quantified by A^, but we also derive limits on erg, the amplitude of matter fluctuations on 8/i"^Mpc scales. 
A higher value for t produces more overall damping of the CMB temperature signal, making it somewhat 
degenerate with the amplitude, A^, and therefore as- The value of erg also affects the height of the acoustic 
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Fig. 3. — Constraints on ACDM parameters from the five-year WMAP data. The two-dimensional 68% 
and 95% marginalized limits are shown in blue. They are consistent with the three- year constraints (grey). 
Tighter limits on the amplitude of matter fluctuations, cts, and the cold dark matter density ^Ich^, arise from 
a better measurement of the third temperature (TT) acoustic peak. The improved measurement of the EE 
spectrum provides a 5(t detection of the optical depth to reionization, r, which is now almost uncorrelated 
with the spectral index Us- 



peaks at small scales, so information is gained from both temperature and polarization. The five-year data 
give (78 = 0.796 ± 0.036, slightly higher than the three-year result, driven by the increase in the amplitude 
of the power spectrum near the third peak. The value is now remarkably consistent with new measurements 
from weak lensing surveys, as discussed in Section [3.21 



3.1.1. Reionization 



Our observations of the acoustic peaks in the TT and TE spectrum imply that most of the ions and 
electrons in the universe combined to make neutral hydrogen and h elium at z c± l l OO. O bservations of 
quasar spectra show diminishing Gunn-Peterson troughs at z < 5.8 |Fan et ahlboool l200lh implying that 
the universe was nearly fully ionized hy z — 5.7. How did the universe make the transition from being nearly 
fully neutral to fully ionized? The astrophysi cs of reionization has been a very active area of research in the 
past decade. Several recent reviews (Barkana fc Loebl 20061 : Fan et al. 2006t Furlanetto et al. 2006t Meiksin 



20071 ) summarize the current observations and theoretical models. Here, we highlight a few of the important 
issues and discuss some of the implications of the WMAP measurements of optical depth. 

What objects reionized the universe? While high redshift galaxies are usually considered the most likely 
source of reionization, AGNs may also have played an important role. As galaxy surveys push towards ever 
higher redshift, it is unclear whether the known population of star forming galaxies at z ~ 6 could have 
i onized the universe (see e.g.. lBunker et al.l (|2007l l). The EE signal clearly seen in the WMAP five-year data 
( Nolta et al. ( 2008h : §2) impHes an optical depth, r ~ 0.09. This large optical d epth suggests tha t higher 
redshift galaxies, perhaps the low luminosity sources appearing in z > 7 surveys (|Stark et al.ll2007r ). played 
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an important r ole in reionization. While the known population of A GNs can not be a significant source 
of reionization jB olton fc Haehnelt 2007 : Srbinovskv fc Wvith e''2007'). an early generation of supermassive 
black holes could have played a role in reionization (jRicotti fc Ostriker„2004ii ; iRicotti et al..i2007l ). This early 
reionization would also have an impact on the CMB. 
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Fig. 4. — Constraints from the five-year WMAP data on ACDM parameters (blue), showing marginahzed 
one-dimensional distribution s and two-dimension al 68% and 95% limits. Parameters are consistent with the 
three-year limits (grey) from lSpergel et al.l (j2007[ ). and are now better constrained. 
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Fig. 5. — Left: Marginalized probability distribution for Zreion in the standard model with instantaneous 
reionization. Sudden reionization at z = 6 is ruled out at 3.5tT, suggesting that reionization was a gradual 
process. Right: In a model with two steps of reionization (with ionization fraction x^, at redshift z^, followed 
by full ionization at z = 7), the WMAP data are consistent with an extended reionization process. 



Mos t of our observational constraints probe the end of the epoch of reionizatio n. Observations of z > 6 



2001 



2001 



Fan et allliooi : Iwillott et allboOTh find that the Lyman-a 
optical depth rises rapidly. Measurements of the afterglow spectrum of a gamma ray burst at z = 6.3 



quasars (jBeckei 



Diorgovski et al 



Totani et al 



20061) suggest that universe was inqstly ionized at z = 6. 3 Lyman alph a emitter surveys 



(jXaniguchi et al.l boOSt [Malhotra k Rhoadsl I2OO6I : IXashikawa et al.l l2006t live et al.l [iooel : lota et all lioOTh 
imply a significant ionized fraction at z = 6.5. The int erpretation that there is a sudden cha nge in the 



properties of the IGM remains a subject of active debate (jBecker et al 



2007 



Wvithe et al.ll2008r ) 



The WMAP data place new constraints on the reionization history of the universe. The WMAP data best 
constrains the optical depth due to reionization at moderate redshift (z < 25) and only indirectly constrains 
the redshift of reionization. If reionization is sudden, then the WMAP data implies that Zroion = 11.0 ± 1.4, 
shown in Figure[5l and now excludes Zreion = 6 at more than 99.9% CL. The combination of the WMAP data 
implying that the universe was mostly reionized at z ~ 11 and the measurements of rapidly rising optical 
depth at z ~ 6 — 6.5 suggest that reionization was an extended process rather than a sudden transition. Many 
early studies of reionization envisioned a rapid transition from a neutral to a fully ionized universe occurring 
as ionized bubbles percolate and overlap. As Figure [5] shows, the WMAP data suggests a more gradual 
process with reionization beginning perhaps as early as z ~ 20 and strongly favoring z > 6. This suggests 
that the universe underwent an extended period of partial reionization. The limits were found by modifying 
the ionization history in CAMB to include two steps in the ionization fraction at late times (z < 30) : 



the first at Z£_with ionization fraction cCr 



Chiu et al 



2006 



the second at z = 7 with x^. — 1. Several studies ( Cen 



2003; 



2003t IWvithe fc Loeblbooil: iHaiman fc Holder] [2OO3I : IVoshida et allbool Ichoudhurv fc Ferrara 



^ Iliev et al. 20071 : Wvithe et aL[2008| ) suggest that feedback produces a prolonged or perhaps even 



multi-epoch reionization history. 



While the current WMAP data constrain the optical depth of the universe, the EE data does not 
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Fig. 6. — Effect of foreground treatment and likelihood details on ACDM parameters. Left: The number 
of bands used in the template cleaning (denoted 'T'), affects the precision to which r is determined, with 
the standard KaQV compared to QV and KaQVW, but has little effect on other cosmological parameters. 
Using maps cleaned by Gibbs sampling (KKaQV (G)) also gives consistent results. Right: Lowering the 
residual point source contribution ('lower ptsrc') and removing the marginalization over an SZ contribution 
('No SZ') affects parameters by < 0.4(7. Using a larger mask ('80% mask') has a greater effect, increasing 
flbh^ by 0.5(7, but is consistent with the effects of noise. 



yet provide a detailed constraint on the reionization history. With more data from WMAP and upcoming 
data from Planck, the EE spectrum will b e gin to place stronger co nstraints on the details of reionization 



(jKaplinghat et al 



2003; 



Holder et al 



200S 



Mortonson Sz Hu 



200: 



3). These measurements will be supple- 



mented by me asurements of the Ostriker-Vishniac effect by high resolution CMB experiment s which is sensi- 
tive to J njdt (Ostriker fc Vishniac 1986t Jaffe &: Kamionkowski 1998t Gruzinov fc Hu IQQSh . and discussed 



m e.g., 



Zhang et all (|2004f ). 



3.1.2. Sensitivity to foreground cleaning 



As the E-mode signal is probed with higher accuracy, it becomes increasingly important to test how 
much the constraint on r, Zroion, and t he other co s molog ical parameters, depend on details of the Galactic 
foreground removal. Tests were done in lPage et al. ( 2007 ) to show that r was insensitive to a set of variations 
in the dust template used to clean the maps. In Figure[6]we show the effect on ACDM parameters of changing 
the number of bands used in the template-cleaning method: discarding Ka band in the 'QV combination, or 
adding W band in the 'KaQVW combination. We find that r (and therefore Zrcion) is sensitive to the maps. 
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but the dispersion is consistent with noise. As expected, the error bars are broadened for the QV combined 
data, and the mean value is r = 0.080±0.020. When W band is included, the mean value is t = 0.100±0.015. 
We choose not to use the W band map in our main analysis however, because there appears to be excess 
po wer in the cleaned m ap at £ = 7. This indicates a potential systematic error, and is discussed further 
(j2008l ). The other cosmological parameters are only mildly sensitive to the number of 



m 



Hinshaw et al 



bands used. This hiK hlights the fact that r is no longer as strongly correlated with other parameters, as in 
earlier WMAP data (jSpergel et al.n2003t 120071 ). notably with the spectral index of primordial fluctuations, 
Us (Figure [3]). 



We also test the parameters obtained using the Gibbs-cleaned maps described in Dunklev et al. ( 2008h 
and Section r2. 1.11 which use the K, Ka, Q, and V band maps. Their distributions are also shown in Figure 
m and have mean r — 0.100 ± 0.018. This is less than la higher than the KaQV template-cleaned maps but 
uses an independent method. The other cosmological parameters are changed by less than 0.3cr compared to 
the template-cleaned results. This consistency gives us confidence that the parameter constraints are little 
affected by foreground uncertainty. The difference in central values from the two methods indicates an error 
due to foreground removal uncertainty on r, in addition to the statistical error, of 0.01. 



3.1.3. Sensitivity to likelihood details 



The likelihood code used for cosmological analysis has a number of variable components that have 
been fixed using our best estimates. Here we consider the effect of these choices on the five-year ACDM 
parameters. The first t wo are the t r eatme nt of the residual point sources, and the treatment of the beam 
error, both discussed in lNolta et al.l (j2008l ). The multi- frequency data are used to estimate a residual point 



source amplitude of Aps — 0.011 ± 0.001 /iK^sr, which scales th e expected cont r ibution to the cros s 



powe r 



spectra of sources below our detection threshold. It is defined in lHinshaw et al.l (|2007( ): iNolta et al.l (j2008r ). 
and is marginalized over in the likelihood code. The estimate co mes from QVW data, w hereas the VW data 
give 0.007 ± 0.003 /xK^sr, both using the KQ85 mask described in lHinshaw et al.l (j2008l ). The right panels in 
Figure [6] shows the effect on a subset of parameters of lowering Aps to the VW value, which leads to a slightly 
higher n^, Jlc/i^ and its, all within 0.4cr of the fiducial values, and consistent with more of the observed high-^ 
signal being due to CMB rather than unresolved point sources. We also use Aps = 0.011 /iK^sr with no 
point source error, and find a negligible effect on parameters (< O.lcr). The beam window function error is 
q uantified by ten modes , and in the standard treatment we marginalize over them, following the prescription 
in Hinsha w et al.l (j2007l ). We find that remov ing the beam error also has a negligible effect on parameters. 
This is discussed further in lNolta et al. ( 20081 ) . who considers alternative treatments of the beam and point 



source errors. 



The next issue is the treatment of a possible contribution from Sunyaev Zeldovich fluctuations. We 
account for the SZ effect in the same way as in the three-year a nalysis, marginalizing ov er the amplitude 
of the contribution parameterized by the Komatsu-Seljak model ( Komatsu &: Seljak 2002h . The parameter 
Asz is unconstrained by the WMAP data, but is not strongly degenerate with any other parameters. In 
Figure [6] we show the effect on parameters of setting the SZ contribution to zero. Similar to the effect of 
changing the point source contribution, the parameters depending on the third peak are slightly affected, 
with a < 0.25(7 increase in rig, flch^, (Ts and similar decrease in baryon density. 

Another choice is the area of s ky used for cosmological interpretation, or how much we mask out to 
account for Galactic contamination. iGold et al.l (|2008i) discuss the new masks used for the five-year analysis. 
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with the KQ85 mask used as standard. We test the effect of using the more conservative KQ80 mask, and 
find a more noticeable shift. The quantity fli,h^ is increased by 0.5a, and Us, flch'^ and ag all decreased 
by ~ 0.4fT. This raised concern s that the KQ85 mask contains residual foreground contamination, but as 
discussed in iNolta et al.l (|2008r ). this shift is found to be consistent with the effects of noise, tested with 
simulations. We also confirm that the effect on parameters is even less for ACDM models using WMAP with 
external data, and that the choice of mask has only a small effect on the tensor amplitude, raising the 95% 
confidence level by ~5%. 

Finally, we test the effect on parameters of varying aspects of the low-^ TT treatment. These are 
discussed in Appendix B, and in summary we find the same parameter results for the pixel-based likelihood 
code compared to the Gibbs code, when both use £ < 32. Changing the mask at low-i to KQ80, or using 
the Gibbs code up to € < 51, instead of ^ < 32, has a negligible effect on parameters. 



3.2. Consistency of the ACDM model with other data sets 

While the WMAP data alone place strong constraints on cosmological parameters, there has been 
a wealth of results from other cosmological observations in the last few years. These observations can 
generally be used either to show consistency of the simple ACDM model parameters, or to constrain more 
complicated models. In this section we compare a broad range of current astronomical data to the WMAP 
A CDM model. A s ubset of the data is used to place combined constraints on extended cosmological models 



Komatsu et al.l (j2008l ). For this subset, we describe the methodology used to compute the likelihood for 



each case, but do not report on the joint constraints in this paper. 



3.2.1. Small-scale CMB measurements 



A number of recent CMB experiments have probed smaller angular scales than WMAP can reach and 
are therefore more sen s itive t o the higher order a coustic oscillations and the details of recombination (e.g., 
Hu fc Sugivamal |l994l . Il996[) : IHu fc Whit3 |l997t )V Since the three-year WMAP analysis, there have been 
new temperature results from t he Arcminute Cosmol ogy Bolometer Array Receiver (ACBAR), both in 2007 
( Kuo et al. 2007) and in 2008 ( Reichardt et al. 2008h . They have measured the angular power spectrum at 
145 GHz to 5' resolution, over ^--^600 deg^. Their results are consistent with the model predicted by the 
WMAP five-year data, shown in Figure [3 although ACBAR is calibrated using WMAP, so the data are not 
completely independent. 

Figure [7] also shows data from the BOOMERANG, CBI and VSA experiments, which agree well with 
WMAP. There have also been new observatio ns of the CMB polar ization from two ground-based experiments, 
QUaD, operating at 100 GHz and 150 GHz |P. Ade et al.ll2007h . and CAPMAP, at 40 GHz and 1 GH z 



(jBischoff et al.l 120081 ). Their measurements of t he EE po wer spectru m are shown in 



together with detect ions already made since 2005 ( iLeitch e t al. 2005; Si evers et al 



Montrov et al. 2006h . and are all consistent with the ACDM model parameters 



2007 



Noltaetal 



( 20081). 



Barkats et aPboOSl : 



In our combined analysis in lKomatsu et al.l (|2008f ) we use t wo different dat a combinations. For the first 
we combine four data sets. This includes the 2007 ACBAR data ( Kuo et al ]l2007h . using 10 bandpowers in the 
range 900<£<2000. The values a nd errors were o btain ed from the ACBAR web site. We also incl ude the three 
external CMB data sets used in ISpergel et al.l ( 2003 ): the Cosmic Background Imager (CBI; Mason et al 
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Fig. 7. — The best-fit temperature angular power spectrum from WMAP alone (red), is consistent with data 
from recent small-scale CMB experiments: ACBAR, CBI, VSA and BOOMERANG. 



(2003) 



Sievers et al 



(|2003h : IPearson et al.l t003i\: 



Readhead et al 



J 2004hl. the Very Small Array (VSA ; 



Dickinson et all |2004l )) and BOOMERANG |Ruhl et al.lboosl : iMontroy et al]|2006l : Ipiacentini et al.ll2006[ l. 



As in the three-year release we only use bandpowers that do not oyerlap with the signal-dominated WMAP 
data, due to non-triyial cross-correlations, so we use seyen bandpowers for CBI (in the range 948 < £ < 1739), 
five for VSA (894 < £ < 1407) and seven for BOOMERANG (924 < £ < 1370), using the lognormal form 
of the likelihood. Constraints are also found by combining WMAP with the 2008 ACBAR data, using 16 
bandpowers in the range 900 < £ < 2000. In this case the other CMB experiments are not included. We do 
not use additional polarization results for parameter constraints as they do not yet improve limits beyond 
WMAP alone. 



3.2.2. Baryon Acoustic Oscillations 



The acoustic pea,k in the galaxy correlation function is a prediction of the adiabatic cosHiologic al model 
|Peebles fc Yulll970l : [Sunvaev fc Zel'dovichl [l97ol : iBond fc Efstathioul llisi IHu fc Sugiyamalll996h . It was 
first detected using t he SPSS luminous red g alaxy (LRG) survey, using the brightest class of galaxies at mean 
redshift z = 0.35 bv lEisenstein et al.l ([20051). The peak was detected at 100/i~^Mpc separation, providing a 
standard ruler to measure the rati o of distances to z = .35 and the CMB at z — 1089, and the absolute 
distance to z — 0.35. More recently [Percival et all (|2007b[ l have obtained a stronger detection from over half 
a million SDSS main galaxies and LRGs in the DR5 sample. They detect baryon acoustic oscil lations (BAO) 



with ov er 99% confidence. A combined analysis was then undertaken of SDSS and 2dFGRS bv lPercival et al 
(j2007a[ ). They find evidence for BAO in three catalogs: at mean redshift z — 0.2 in the SDSS DR5 main 
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Fig. 8. — Baryon Acoustic Oscillations expected for the best-fit ACDM model (red lines), compared to 
BAO in galaxy power sp ectra calculated from (left) combined SDSS and 2dFGRS main galaxies, and (right) 
SDSS LRG galaxies, bv IPercival et al.l (|2007a|). The observed and model p ower spectra have been divided 
by P(fc)smooth, a smooth cubic sphne fit described in iPercival et alj (|2007a[) . 



galaxies plus the 2dFGRS galaxies, aX z = 0.35 in the SDSS LRGs, and in the combined catalog. Their data 
are shown in FigurelU together with the WMAP best- fit model. The BAO are sh own by dividing the ob served 
and model power spectra by P(fc)smooth, a smooth cubic spline fit described in IPercival et al.l (j2007al ). The 
observed power spectra are model-dependent, but were calculated using f2„i = 0.25 and h = 0.72, which 
agrees with our maximum-likelihood model. 

The scale of the BAO is analyzed to estimate the geometrical distance measure at z = 0.2 and z — 0.35, 

Dv{z)^[{l + zfD\cz/H{z)Y/\ (9) 

where Da is the angular diameter distance and H{z) the Hubble parameter. They find rg/ Dv{0.2) = 0.1980± 
0.0058 and rs/Z)y(0.35) = 0.1094 ± 0.0033. Here is the comoving sound horizon scale at recombination. 
Our ACDM model, using the WMAP data alone, gives r ^15^(0.2) = 0.1946 ± 0.0079 and r ^15^(0.35) = 
0.1165 ± 0.0042, showing the consistency between the CMB measurement at z = 1089 and the late-time 
galaxy clustering. However, while the z = 0.2 measures agree to within la, the z = 0.35 measurements have 
mean values almost 2cr apart. The BAO constraints are tighter then the WMAP predictions, which shows 
that they can improve upon the WMAP parameter determinations, in particular on JIa and flmh^- 



described in 



Percival et al 



In iKomatsu et al. ( 2008) the combined bounds from both surveys are used to constrain models as 



( 2007a[) . adding a likelihood term given by — 21nL ~ X'^C ^X, with 
X^ ^ [r,/Dv{0.2) - 0.1980, r,/L>y(0.35) - 0.1094] 



(10) 



and 



35059^12 



24031, C — 108300, including the correlation between the tw o measurements. 



Komatsu et al.l (|2008l ) also consider constraints using the SDSS LRG limits derived by lEisenstein et al 
|2005h . using the combination 

' (11) 



for z 



A{z) = Dv{z)^nr^Hl/cz 
0.35 and computing a Gaussian likelihood 

- 21nL = (A - 0.469(n,/0.98)-°-^^)V0.017^. 



(12) 
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3.2.3. Galaxy power spectra 



We can compare the predicted fluctuations from the CMB t o the shape of galaxy power spectra, in 
addition to the scale of acoustic o scillations fe.g. jEisenstein fc Hul (liggsl)). The SPSS galaxy power spectum 
from DR3 |Tegmark et al.ll2004l ) and the 2dFGRS spectrum ([ColeeLalJ l2005|) were shown to be in gooc 
agreement with the WMAP three- year data, and used to place tighter constraints on cosmological models 
(jSpergel et al.ll2007n . but there was some tension between the preferred values of the matter density {Vl„i ~ 
0.236 ± 0.020 for WMAP combined with 2dFGRS and 0.265 ± 0.030 with SDSS). Recent studies used 
photometric redshifts to estimate the galaxy power spectrum of LRGs in the range 0.2 < z < -6 fro m the 
SDSS fourth data release (DR4), finding »^ = 0.30 ± 0.03 (for h = 70, IPadmanabhan et all |2007h ) and 
n^h = 0.195 ± 0.023 for h = 0.75 (iBlake et al.lboOTi ). 



More prec i se me asurements of the LRG power spectrum were obtained from rcdshift measurements: 
Tegmark et"al] |2006h used LRGs from SDSS DR4 in the range O.Ol/i/Mpc < A: < 0.2ft./Mpc combined 
with the three-year WMAP data to place strong constraints on cosmological models. However, there is a 
disagree ment between the matt er density predicted using different minimum scales, if the non-linear modeling 
used in iTegmark et al.l (j2006l ) is adopted. Using the three-year WMAP data combined with the LRG 
spectrum we find Vim = 0.228 ± 0.019, using scales with k < O.lh Mpc-\ and = 0.248 ± 0.018 for 
k < 0.2h Mpc~^. Th ese constraints are ob tained for the 6 parameter ACDM model, following the non- 
linear prescription in iTegmark et al.ll2006h . The predicted galaxy power spectrum Pg{k) is constructed 
from the 'dewiggled' linear matter power spectrum Pm{k) using Pg{k) = 6^Pm(fc)(l -I- Qk'^)/{1 -f 1.4A:), and 
the parameters b and Q are marginalized over. The de wiggled spectrum is a weighted sum of the matter 
power spectrum at z = and a smooth power spectrum ( Tegmark et al. 2006h . Without this dewiggling, we 
find Jim — 0.246 ± 0.018 for k < 0.2h Mpc~^, so its effect is small. We also explored an alternative form 
for the bias, motivated by third-orde r perturbation theory analysis, with Pg{k) = b'^P^l{k) + N (see e.g.. 



200lh:IS chulz fc Whitel I2OO6I )). Here the non-linear matter power spectrum using the Halofit 



Seljak 
model 

Marginalizing over h and N we still find a discrepancy with scale, with = 0.230 ±0.021 using scales with 
k < O.lh Mp c~^, and 0^ = 0.249 z t 0.018 for k < 0.2h Mpc^^. Constraints using this bias model are also 



Smith et al.ii2003l ) and iV is a constant accounting for non-linear evolution and scale-dependent bias. 



considered in^Hamann etal 



( 2008a). These fin dings are consistent with results obtained from the DR5 main 



galaxy and LRG sample (jPercival et al.ll2007c[ ). who argue that this shows evidence for scale-dependent bias 
on large-scales, which could explain the observed differences in the early SDSS and 2dFGRS results. While 
this will likely be resolved with future modeling and obser vations, we choose not to use the galaxy power 
spectra to place joint constraints on the models reported in Komatsu et al. ( 2008[) . 



3.2.4. Type la Supernovae 



In the last decade Type la supernovae have become an important cosmological probe, and have provided 
the first direct evidence for the acceleration of the universe by measuring the luminosity distance a s a function 
of red s hift. The observed di n iness of high redshift su pernovae {z ~ 0.5) was first measured bv lRiess et al 



(1998) 



Nobili et al 



Schmidt et al 
2OO5I ): 



(1998) 



Perlmutter et al 



Krisciunas ^ 



j 1999 ), confirmed wit h more recent measu rements including 



et all |2005[ l: lciocchiatti et al.l |2006[ ): Ustier et all |2006h . and extended to 
higher redshift by iRiess et al.l ( 20041 ) who found evidence for the earlier deceleration of the universe. The 
sample of high redshift supernovae has grown by over 8 since the t hree-y ear WMAP analysis. Recent 



HST measurements of 21 new high redshift supernova bv iRiess et al.1 (|2007| ) include 13 at z > 1, allowing 
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Fig. 9. — The red line shows the luminosity distance relationship predicted for the best-fit WMAP- 
only model (right column in Tab le 2). The points show binned supernova observations from the Union 



compilation jKqwalski et ahl boosh. including high red shift SNIa from HST |Riess et ahlboOTl ). ESSENCE 



(jMiknaitis et al.l 120071 ). and SNLS (jAstier et al.ll2006( ). The plot shows the deviation of the luminosity 



distances from the empty universe model. 



the measurement of the Hubble expansion H{z) at distinct epochs and strengthening the evidence for a 
period of deceleration followed by acceleration. The ESSENCE Supernova Survey has also recently reported 
results from 102 supernovae discove red from 2002 to 200 5 using the 4-m Blanco Telescope at the Cerro 
Tololo Inter-American O bservatory (jMiknaitis et al.l 120071 ). of which 60 are used for cosmological analysis 
(jWood-Vasev et al.ll2007|). A comb i ned c osmological analysis was performed of a subset of the c omplete 
supernova data set bv iDavis et al.l (|2007l ) using the MCLS2k2 light curve fitter (jjha et al.ll2007l) . More 
recently, a b roader 'Union' compil ation of the currently observed SNIa . including a new nearby sample, was 
analysed bv lKowalski et all |2008l ) using the SALTl hght curve fitter |Guv et al.lboosh . 



We confirm in Figure [9] that the recently observed supernovae are consistent with the ACDM model, 
which pr edicts the luminosity distance fith as a function of redshift and is compared to the Un ion combined 
data set ( Kowalski et al. 2008h . In the cosmological analysis described in Komatsu et al. I liooa . the Union 
data is use d, consisting of 307 supernovae that pass various selection cuts. These include supernovae observed 
using HST lRiess et al.l ^OOtI ) , from the ESSENCE survey jMiknaitis et al.ll2007t ). and the Supernova Legacy 
Survey (SNLS, Astier et al. ( 2006h ). For each supernova the luminosity distance predicted from theory is 
compared to the observed value. This is derived from measurements of the apparent magnitude m and 
the inferred absolute magnitude Af, to estimate a luminosity distance ^obs = 51og[(ii(z)/Mpc] + 25. The 
likelihood is given by 



21nL = 2_^[^j.obsA{zi) - f^thA^i,Mo)] /fJohs^i 



(13) 



summed over all supernovae, where a single absolute magnitude is marginalized over ( Lewis fc Bridlel 2002 ) 
and (Jofcs is the observational error accounting for extinction, intrinsic redshift dispersion, K-correction and 
light curve stretch factors. 
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Data 


Parameter 


Lensing limits 


5-year WMAP limits 


CFHTLS Wide 

100 Sq Deg 
COSMOS 2D 
COSMOS 3D 


a8(a„/0.25)0-64 
a8(a„/0.24)0-^9 
a8(r!™/0.3)"-4« 
f78(17™/0.3)0-44 


0.785 ±0.043 
0.84 ±0.07 
0.81 ±0.17 


0.814 ±0.090 
0.832 ±0.088 
0.741 ± 0.069 
0.745 ±0.067 



Table 3: Measurements of combinations of the matter density, flrn, and amplitude of matter fluctuations, as, from 
weak lensing observations ( Fu et al. 20081 : Benjamin et al. 2007 : Massev et al. 2007 ). compared to WMAP. 



3.2.5. Hubble constant measurements 



The WMAP estimated value of the Hubble constant, Ho — 71.9^2 7 km s ^ Mp c ^ , assuming a flat ge - 
ometry, is consistent with the HST Key Project measurement of 72±8 km s^^ Mp c~^ (Freedman et al I2OO1I) . 
It also agrees within Icr with measurements from gravitationally lensed systems ( Koopmans et al.l 2003), S Z 



and X-ray observations ( Bonamente et al. 2006h . Cepheid distances to nearby galaxies faiess e t al. 



the d istance to the Maser-host galaxy NGC4258 as a calib rator for the Cepheid distance scale (jMacri et al 



2005) 



2006 ). and a new measure of the TuUy-Fisher zero-point ( Masters et al. 20061 ). the latter two both giving 
Hq = 74±7 km s~^ Mpc~^. L ower measures are fav ored by a compilation of Cepheid distance measurements 
for ten galaxies using HST bv lSandage et al.l (|2006r ) (i?o = 62 ± 6 km s~^ Mpc~^), and new measurements 
of an eclipsing binar y in M33 which reduce the Key Project measurement to Hq — 61 km s~^ Mpc 



( Bonanos et al. 20061 ). Higher measures are found using revised par allaxes for Cepheids ( van Leeuwen et al 



200 



3), raising the Key Project value to 76 ± 8 km s"^ Mpc"^ In iKomatsu et al.l J2008h the Hubble Con 



stant measurements are included for a limited set of parameter constraints, using the lFreedman et al.l (|200ll ) 
measurement as a Gaussian prior on Hq. 



3.2.6. Weak Lensing 



Weak gravitati onal lensing is pr oduced by the distortion of galaxy images by the mass distribution along 
the line of sight (see Refregier ( 2003h for a review). There have been sig nificant advan c es in i ts measurement 
in recent years, a nd in the underst anding of systematic effects (e.g., iMassev et al.l (|2007| )). and intrinsic 



alignment effects (iHirata et al. 2007), making it a valuable cosmologica l probe complementary to t he CMB. 
Early results by the PCS (jHoekstra et al.ll2002l ). VIRMOS- DESCART Jvan Waerbeke et al.l l2005h. and the 



Canada-France-Hawaii Telescope Legacy Survey fCFHTLS |Hoekstra et al.l (|2006r i) lensing surveys favored 



higher amplitudes of mass fluctuations than preferred by WMAP. H owever, new m easurements of the two- 
point correlation functions from the third year CFHTLS Wide survey (|Fu et al l2008h favor a lower amplitude 
consistent with the WMAP measurements, as shown in Table |3l This is due to an improved estimate of 



the galaxy r edshift distribution from CFHTLS-Deep (|Ilbert et al.ll2006l ). compared to that obtained from 



photometric redshifts from the small Hubble Deep Field, which were d ominated by systemat ic errors. Their 
measured signal agrees with results from the 100 Square Degree Survey ([Benjamin et al ] l2007[l . a compilation 
of data from the earlier CFHTLS - Wide, RCS and VIRMOS-DESCAT surveys, together with the GABoDS 
survey (iHetterscheidt et al.l 120071 ). with average source redshift z ~ 0.8. Both these analyses rely on a 
two-dimensional measurement of the shear field. Cosmic shear has also been measured in two and three 



dimensions by the HST COSMOS survey (jMassev et al.l 120071 ). using redshift information to providing an 
improved measure of the mass fluctuation. Their measures are somewhat higher than the WMAP value, as 
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shown in Table [3l although not inc onsistent. Weak lensing is also produced by the d istortion of the CMB by 
the intervening mass distribution (IZaldarria ga &: Seliak 1999 : Hu fc Okamotol 2002 ). and can be probed by 
measuring the correlation of the lensed CMB w ith tracers of large scale structure. Tw o recent analyses have 
found evidence for the cross-correlated signal ( Smith et al. 2007 : Hirata et al. 2008h . both consistent with 
the five-year WMAP ACDM model. They find a 3.4g d etection of the correlation between the three-year 
WMAP data and NVSS radio sources (jSmith et al.ll2007l ). and a corr elation at the 2.1g level of significance 
between WMAP and data from NVSS and SDSS LRGs and quasars |Hirata et allboosl ). 



3. 2. 7. Integrated Sachs- Wolfe Effect 



Correlation between large-scale CMB temperature fluctuations and large-scale structure is expected in 
the ACDM mode l due to the ch ange in gravitational potential as a function of time, and so provides a test 
for d ark energy (iBouehn et al. Evidence of a c orrelation was found in the first-year WMAP data 



(e.g., iBoughn fc CrittendenI (|2004f ): 



Nolta et al. 



(|2004l )). Two recent analyses combine recent large-scale 
structure data (2MA SS, SPSS LRG s, SDSS q uasars and NVSS radio so urces) with the WMAP three-year 
data, finding a 3.7a ( Ho et al. 2008h and 4ct ( Giannantonio et al. 2008h detection of ISW at the expected 
level. Other re cent studies using individual data sets find a correlation at the level expected with the SDSS 
DR4 galaxies JCabre et allbood). at high redshift with SDSS qusar s |Giannantonio et al.ll200^ . and with 



the NVSS radio galaxies |Pietrobon et aLlliool iMcEwen et al.ll2007 l ) . 



3.2.8. Ly-a Forest 



The Lya forest seen in quasar spectra probes the underlying matter distribution on small scales (j Ranch 

19981 ). However, the relationship between absorption line structure and mass fluctuations must be fully 
understood to be used in a cosmological analysis. The power spectrum of the Lyg forest has been used 



to constrain the shape and amplitude of the primo rdial power spectrum (jViel et al.l 120041 : iMcDonald et al 



2005t ISeliak et al.ll2005l : iDesiacaues &: Nusseij|2005l ) , and recent results combine the th ree-year WMAP data 
with the power spec trum obtained fro m the LUQAS sample of VLT-UVES spectra ( Viel et al. 2006[ ) and 



SDSS QSO spectra |Seliak et al.ll200fih . Both grou ps found results s uggesting a higher value for erg than 



consistent with WMAP. However, measurements bv iKim et al.l (j2007h of the probability distribution of the 
Lyg flux have been compared to simulations with different cosmological parameters and thermal histories 
(jBolton et al.ll2007l ) . They imply that the temperature-density relation for the IGM may be close to isother- 
mal or inverted, which would result in a smaller amplitude for the power spectrum than previo usly inferred, 



more i n line with the flve-year WMAP value of as = 0.796±0.036. Simulations in larger boxes bv lTvtler et al 



([20071) fail to match the distribution of flux in observed spectra, providing further evidence of disagreement 
between simulation and observation. Gi v en th ese uncertainties, the Lya forest data are not used for the 
main results presented in lKomatsu et aL ( 2008). However , constraints on the running of the spectral index 
are discussed, using data described in ISeliak et al.l \200(h . With more data and further analyses, the Lya 
forest measurements can potentially place powerful constraints on the neutrino mass and a running spectral 
index. 
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3.2.9. Big-Bang Nucleosynthesis 



WMAP measures the baryon abundance at decoupling, with ilt,h^ = 0.02273 ± 0.00062, giving a baryon 
to photon ratio of 7yio(WMAP) — 6.225 ± 0.170. Element abundances of deuteri um, helium, and lithium 
also depend on the baryon abundance in the first few minutes after the Big Bang. ISteigmanI (j2007l ) reviews 
the current status of BBN measurements. Deuterium measurements provide the strongest test, a nd ar e 
consistent with WMAP, giving 7710 (-D) — 6.0 ± 0.4 based on new measurements bv lO'Meara et al.l (j2006l ). 
The '^He abundanc e is more poorly constrained at ?7io(^He) — b.G'^l'l from the measure of 1/3 — 1.1 ± 0.2 by 
Bania et all |2002h . The neutral lithium abundance, measured in low-metallicity stars, is two times smaller 
than the CMB prediction, ?7( Li) — from measures of the l ogarithmic abundance, [Lijp = 12 + logmfLi/H) , 



is the range [Li]p ^ 2.1 — 2.4 ( Charbonnel fc Primasl2005 : Melendez fc Ramire: 



Bird et al 



2004 



Asplund et al.ll2005l). 



Thes e measurements c ould be a signature of new early universe physics, e.g., (jCoc et al.ll2004 : 
2005t 



Richard et al 



Jedamzik 2004 ), with r ecent at tempts to simultan e ously fi t both the "^Li an d ^Li abunda n ces b y 



|2007h : IjedamzikI |2008ah : ICumberbatch et al.l |2007l ): IjedamzikI |2008b[) : Ijittoh et all (|2008V 



The discrepancy could also be due to systemati cs, destruction of l ithium in an earlier generation of st ars, 
or uncertainties in the stellar temparature scale ( Fields et al. 20051 : Steigma ^ l2006HAsDlund et al.ll200,5h. A 



possible sol ution has been proposed using observations of stars in the globular cluster NGC 6397 (jKorn et al 
2006], |2007|). They find evidence that as the stars age and evolve towards hotter surface temperatures, the 
surface abundance of lithium in their atmospheres drops due to atomic diffusion. They infer an initial 
lithium content of the stars [Li]p = 2.54 ± 0.1, giving 77io(Li) — 5.4 ± 0.6 in good agreement with BBN 
predictions. However, further observations are needed to determine whether this scenario explains the 
observed uniform depletion of primordial ^Li as a function of metallicity. The measured abundance of "^He 
is also low er than predicted , with r jioj^He) = 2.7ln g from a measure of Yp = 0.240 ± 0. 006, incorporatiri g 
data from llzotov fc Thu anI (fiooi); loiive & SkillmaJ ('2004'); ' Cruenwald et al.' hooj) bv ISteigmanI |2007l ). 
However, observations bv .Peimbert et al.. (,20071 ) and lFukugita fc Kawasaki (.2006 ) predict higher values more 
consistent with WMAP. 



3.2.10. Strong Lensing 



The number of strongly lensed quasars has the potenti al to probe cosmology, as a dark ene rgy dominated 



universe predicts a large number of gravitational lenses (| Turner 



1990; 



radio band surve y has a large statistical sample of radio le nses (Mvers et al 
York et al.ll2005l ). yielding estimates for f^A ^ 0.72 - 0.78 |Mitchell et al 



2003 



Fukugita et al.lll990l). The CLASS 



Koopmans 



et al. 



2003; 



2005l : IChaell2007l ). lOguri et al 



( 20081 ) have recently anal yzed the large statistical lens sample from the Sloan Digital Sky Quasar Lens 



-o.o6(syst.). 



Search (jOguri et al.l 120061 ). For a w = -1, flat cosmology, they find VLa = 0.74t°;i[^(stat.)t°-^^' 
These values are all consistent with our best fit cosmology. The abundance of giant arcs also has the 
potential to pr obe the underlying cosmology. However, although recent surve ys have detected larger numbers 



of giant arcs (IGladders et al 



with ACDM (iLi et al 



Hennawi et al 



2007 



2OO1 



I 



2003 



Sand et al 



2005 



Broadhurst h Barkana 



Hilbert et al 



2007a; 



Hennawi et al. 20081 ) than argued to be consistent 



Neto et al. 



20081). numerical simulations (jMeneghetti et al.l 12007 



20071 ) find tha t the lens cross-sections are sensitive to 



the m ass distribution in clusters as well as to the baryon physics (jWambsganss et al 
2007b| ). These effects must be resolved in order to test ACDM consistency. 



2007 



Hilbert et al 
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3.2.11. Galaxy clusters 

Clusters are easily detected and probe the high mass end of the mass distribution, so probe the am- 
plitude of density fluctuations and of large-scale structure. Cluster observations at optical wavelengths 
provided some of the first evide nce for a low density universe with the current preferred cosmological param- 



eters (see e.g.. iFan et alJ (|l997l )). Observers are now using a number of differe nt techniques fo r ident ifying 



cluster samples: large optical samples, X-ray surveys, lensing surveys (see e.g.. IWittman et al.l (j2006l )) and 



Sunyaev-Zeldovich surveys. The ongoing challenge is to determine the selection function and the relation- 



(Lin et al. 20061 Sheldon et al. 2007; Reves et al. 2008; Rvkoff et al 


2008) and the X-rav (.Sheldon et al. 


2001; Reiorich & B6hringei]|2002 ; iKravtsov et al. 2006; 


Arnaud et al. 


2007 


Hoekstra 2007). With large new 


optical cluster samples ( Bahcall et al. 2003t Hsieh et al. 


2005; Miller et al. 


2005; Koester et al. 2007) and X- 


ray samples from ROSAT, XMM-LSS and Chandra (Pierre et al. 20061 Burenin et al. 2007; Vikhlinin et al. 



20081 ). cosmological parameters can be further tested, and most recent results for erg are conv erging on values 
close to the WMAP best fit value of cts = 0.796 i b 0.036. Constrain ts from the RCS survey (jCladders et al 
20071 ) give Vl„. 



O.SOtJ^ ii and gs = 0.70 



+0.27 



level) from S PSS BCG samples. iMantz et al.l (jaOOJfind O^, = 0.27+!^-^g and ctr = 0.7 7 



Rozo et al.l (|2007l ) argues for erg > 0.76 (95% confidence 



7-1-0.06 



+1'% for a flat model 



based on the Jenkins et al. ( 20011 ) mass function and the Reiprich fc Bohringei ( 20021 ) mass- luminosity cali 
bration. Wi th a 30% hig h er zer o-point, Rvkoff et al. ( 20081 ) find that their data are best fit by erg = 0.85 and 



(120071) report ag = 0.92^°J[; for 0„ = 0.24 from their joint CFHTLS/XMM-LSS 



Q.m = 0.24. iBerge et al 
analysis. In iRines et al 
sample of X-ray clusters. For VLm — 0.3, they find cg — 0.84 ± 0.03. 



20071 ) redshift data from SDSS DR4 is used to measure virial masses for a large 



3.2.12. Galaxy peculiar velocities 



With deep large-scale structure surveys, cosmologists have now been able measure j3, the amplitude 
of redshift space distortions as a function of redshift. Combined with a measurement of the bias, 5, this 
yields a determination of the growth rate of structure / = dlnG/dlna = 0h, where G is the gro wth 
factor. For Einstein gravity theories we expect f ~ W with 7 ~ 6/11 fsee iPolarski fc Gannouiil (j2007l ) for 
a more accurate fitting function). Ana l ysis of redshift space distortions in the 2dF galaxy redshift survey 
|Peacock et ahlboOll Fverde et al.|[2002t iHawkins et aLlliooj ) find P = 0.47 ± 0.08 a t z 0.1, consistent 
with the ACDM predictions for the best fit WMAP paramete rs. The Tegmark et al. ( 20061 ) analysis of the 
SDSS LRG sample finds f3 = 0.31 ± 0.04 at z = 0. 35. The IRoss et al.l (l2007l) analvsis of the 2dF-SDSS 



LRG sample find (3 = 0.45 ± 0.05 at z = 0.55. The ICuzzo et al.l |2008l ) analysis of 10,000 galaxies in the 
VIMOS-VLT Deep Survey finds that (3 = 0.70 ± 0.26 at z = 0.8 and infer d\nG/dlna = 0.91 ± 0.36, 
cons istent with the more ra pid growth due to matter domination at this epoch expected in a ACDM model. 



The Ida Angela et al.l (|2008f ) analysis of a QSO sample finds /3 = 0.60_q ]^]^ at z = 1.4 and use the clustering 



+0.09 A 
0.07- ^ 



length to infer the bias. Extrapolating back to z = 0, they find a matter density of flm — 0.25 
second approach is to use objects wit h well-determined distances, such as galaxies and supernoyae, to look 
for deviations from the Hu bble fl ow ( Strauss fc Willick 1995 ; Dekel 200d Zaroubi et al. 2001; Riess et al 



1995 



Haugbolle et"al]l2007h . The IPark fc ParkI (120061) an alvsis of the peculiar velocities of galaxies in the 



SFI sample find crgfij^^ = 0.56^0^21- I Gordon et al.l (j2007l ) correlate pecuhar velocities of nearby supernova 
and find g g = 0.79 zt 0.22. These measureme nts provide an independent consistency check of the ACDM 



model (see Nesseris fc Perivolaropouloa ( 2008 ) for a recent review). 



- 26 - 




Fig. 10. — Two-dimensional marginalized constraints (68% and 95% confidence levels) on inflationary pa- 
rameters r, the tensor-to-scalar ratio, and Ug, the spectral index of fluctuations, defined at ko = 0.002/Mpc. 
One-dimensional 95% upper limits on r are given in the legend. Left: The five-year WMAP data places 
stronger limits on r (shown in blue) than three- year data (grey). This excludes some inflationary models 
including A^"* monomial inflaton models with r ~ 0.27, 0.95 for 60 e-folds of inflation. Right: For 

models with a possible running spectral index, r is now more tightly constrained due to measurements of 
the third acoustic peak. Note: the two-dimensional 95% limits correspond to A(21nL) ~ 6, so the curves 
intersect the r = line at the 2.5a limits of the marginalized distribution. 



4. Extended cosmological models with WMAP 



The WMAP data place tight constraints on the simplest ACDM model parameters. In this section we 
describe to what extent WMAP data constrain extensions to the simple model, in terms of quantifying the 
primordial fluctuation s and determining the composition of the universe beyond the standard components. 



Komatsu et al.l (|2008r ) present constraints for WMAP combined with other data, and offer a more detailed 



cosmological interpretation of the limits. 



4.1. Primordial perturbations 

Tensor fluctuations 



In the ACDM model, primordial scalar fluctuations are adiabatic and Gaussian, and can be described 
by a power law spectrum, 

A^(fc)cx - , (14) 



producing CMB angular power spectra consistent with the data. Limits can also be placed on the amplitude 
of tensor fluctuations, or gravitational waves, that could have been generated at very early times . The y 
leave a distinctive large-sc ale signature in the polarized B-mode of the CMB fe.g.. lBasko fc PolnarevI (|l980( ): 
Bond fc Efstathioul (jl984l )'). that provides a clean way to distinguish them from scalar fluctuations. However, 
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Fig. 11. — Two-dimensional marginalized limits for the spectral index, Us, defined at fco — 0.002/Mpc, and 
the running of the index dris/dhik (marked rirun)- Models with no tensor contribution, and with a tensor 
contribution marginalized over, are shown. In both cases the models are consistent with a power-law spectral 
index, with dUg/dlnk — 0, as expected from the simplest inflationary models. 



we have not yet reached sensitivities to strongly constrain this signal with the polarization data from WMAP. 
Instead we use the tensor contribution to the temperature fluctuations at large scales to constrain the tensor- 
to-s calar ratio r. We define r = A^(fco)/A^(A:o), where is the amplitude of primordial gravitational waves 



(see 



Komatsu et al.l (|2008|)), and choose a pivot scale feg = 0.002/Mpc. 



The WMAP data now constrain r < 0.43 (95% CL). This is an improvement over the three-year limit 
of r < 0.65 (95% CL), and comes from the more accurate measurement of the second and third acoustic 
peaks. The depend ence of the tensor amplitude on the spectral index is shown in Figure 1101 showing the 
Us — r degeneracy (jSpergel et al.ll2007l) : a larger contribution from tensors at large scales can be offset by 
an increased spectral index, and an overall decrease in the amplitude of fluctuations, shown in Table |4l The 
degeneracy is partially broken with a better measure of the TT spectrum. There is a significant improvement 
in the limit on models whose scalar fluctuations can vary with scale, with a power spectrum with a 'running' 
spectral index, 



The limit from WMAP is now r < 0.58 (95% CL), about half the three-year value r < 1.1 (jSpergel et al 
2007h . 



What do these limits tell us about the early universe? For models that predict observable gravitational 
waves, it allows us to exclude more of the parame ter space. The sirnplest inflationary models predict a nearly 
scale-invariant spectrum of gravitational waves (|Grishchuklll97i IStarobinskvlll979h . In a simple classical 
scenario where inflation is driven by the potiential Vi^cj)) of a slowly rolling scalar fleld, the predictions 



- 28 - 



Parameter Tensors Running Tensors+Running 



r 



< 0.43 (95% CL) < 0.58 (95% CL) 

dris/dlnk -0.037 ±0.028 -0.050 ±0.034 

n, 0.986 ±0.022 1.03lt!^:Sj^5 1.087;g;^^3 

0-8 0.7771^1^4? 0.816 ± 0.036 0.800 ±0.041 



Table 4: Selection of cosmological parameter constraints for extensions to the ACDM model including tensors and/or 
a running spectral index. 



( Lvth fc Riottolll99g ) are 



r ~ Aa/N 
I -Us ~ {a + 2)/2N 



(16) 
(17) 



for V{4>) (X 0", where N is the number of e-folds of inflation between the time when horizon scale modes left 
the horizon and the end of inflation. For N=60, the A(/)* model with r ~ 0.27, ~ 0.95 is now excluded 



with more than 95% confidence. An 



model with 



0.13, n, ~ 0.97 is still consistent with the data. 



Komatsu et al.l (|2008l ) discuss in some detail what these measurements, and constraints for combined data 
sets, imply for a large set of possible inflationary models and potentials. 

With r — also fitting the data well, mo dels that do not pre dict an o bservable level of grav itational 
waves, including m ulti-field inflationary models ( Polarski fc Starobin skv 1995 HGarcia-Bellido fc Wands. 1996) . 



D-brane inflation (jBaumann fc McAllisteij 120071 ) , and ekpyrotic or cyclic scenarios (jKhourv et al 
Bovle et al.l 120041 ) . are not excluded if one flts for both tensors and scalars. 



2001 



4- 1-2. Scale dependence of spectral index 



The running of the spectral index has been the subject of some debate in light of WMAP observations, 
with the three- year dat a giving limits of dng Id In k = —0.055 ±0.030, showing some preference for decreasing 
power on small scales (|Spergel et al.l 120071) . Combined with high-^ CBI and VSA CMB data, a negative 
running was preferred at ^ 2a. A running index is not predicted by the simplest inflationary models (see e.g., 
Kosowskv fc Turnej |l995[ )). and the detection of a scale dependence would have interesting consequences 
for early universe models. Deviations from a power law index , and their consequences , have been considere d 
by a number of groups in lig h t of three-year data, i ncluding Easther fc Peiris ( 2006[ ): Kinnev et al. (2006); 
Shafieloo fc Souradeepl (120071 ): I Verde fc PeirisI (l2008l) . using various parameterizations. In this analysis, and 



Komatsu et al.l (|2008( l we consider only a running index parameterized as in Eqn ll5l 



We show in Figure[TT]that with a better determination of the third acoustic peak, coupled with improved 
beam determination, the five-year WMAP data do not significantly prefer a scale-dependent index. The limit 
on the running is dns/dlnk — —0.03 7 ± 0.028 for models with no tensor contribution. The running is anti- 
correlated with the tensor amplitude ( Spergel et al.l[2003 . 2007), so the positive prior on the tensor amplitude 
leads to a more negative running preferred, dus/dlnk = —0.050 ± 0.034, when a tensor contribution is 
marginalized over. Both limits are well within 2a of zero, showing no evidence of departure from a power 
law spectral index. 




Fig. 12. — Marginalized two-dimensional limits (68% and 95%) on the amplitude of possible CDM entropy 
(or isocurvature) fluctuations. The one-dimensional 95% upper limits are given in the legend. Left: Anti- 
correlated fluctuations are tightly constrained, placing limits on curvaton models. Right: Uncorrelated 
fluctuations, corresponding to axion models for dark matter, add less power to the CMB spectrum than the 
anti-correlated case, for a given a, so higher values of ao are allowed (than a_i), by the data. In both cases 
the amplitude is correlated with the spectral index of curvature fluctuations n^, which compensates for the 
large scale power added by the CDM entropy fluctuations. 



4-1.3. Entropy perturbations 



The simplest classical single-field inflation models predict solely adiabatic fluctuations , but entropy (or 

isocurvature) fluc tuations are also pre dicted in a wide range of scenarios, including axions ( Seckel fc Turner 

1985 : lLindelll985h . multi-fleld inflation JPolarski fc Starobinskyll99d:lGarcia-Bemdo fc WandJl99d:lLinde fc Mukhanov 
1997 ) and decay of fiel ds such as the curvaton |Lvth fc WandslboojlMoroi fc Takahashil200ll : lBartolo fc Liddle 
20Q2HLvth et al.ll2003|) . They may be correlated with the adiabatic fluctuations to some d egree, depending on 
the model. Most ph ysical scenarios generate only CDM or baryon entropy fluctuations ([Bond fc Efstathiou 
1987t|PeebQll987l) . with perturbation 

(18) 

Pc ^Pt 



for CDM with density pc. Neutrino modes are also possible (IBucher et al. 2000 ) . It has been known for some 
time that the CMB data cannot be fit by pure entropy fluctuations |stompor et al.lll996l : Langlois fc Riazuelol 
20001 ). but a contribution may be allowed . Several groups have placed limits on a variety of models using 



Keskitalo et al 



the WMAP one- year and three-year data (IPeiris et al.ll2003 : IValiviita fc Muhonenll2003t iBucher et al 



Beltran et al ]l2004l:lDunklev et al.ll200,'Tl: iKurki-Suonio et a 



2005; Lewii 



sl l200fit 



Bean et al.f2006l : iTrotta 



2004; 



2007 



200 



3), finding no stro ng evidence for ent r opy fluctuations. Significant levels hav e been found 



to be consistent with the CMB data |Bucher et al.ll2004 Moodlev et aUboO^ iBean et aLlliooj ). but require 
correlated admixtures of CDM and neutrino isocurvature perturbations that are hard to motivate physically. 
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Parameter 



/3 = -l 



ao 

<^8 



< 0.011 (95% CL) 



0.983 ±0.017 

n 77O+0.039 



< 0.16 (95% CL) 
0.987 ±0.022 
0.777 ±0.038 



Table 5: Subset of cosmological parameter constraints for ACDM models with additional anti-correlated (/3 
or uncorrelated (/3 — 0) entropy fluctuations. 



-1) 



Here , and in Komatsu et al. J 2008 ), we quant ify the relative contributions to the angular power spectrum 
following Beltran et all |2004h : lBean et all |2006h . with 



Ce = {l- a)Cf + aCf + 2l3^Ja{l - a)Cf , (19) 

summing the spectra from curvature fluctuations C^, entropy fluctuations Cf with power spectrum 

A|(fc)53(k - k') = (fc/27r)3(5,(k)5e(k')), (20) 

and a cross-correlation spectrum Cf^ with power spectrum 

A|(fc)(53(k - k') = (fc/2^)3(-7^(k)5,(k')). (21) 

This follows the definition of the curvature perturbation TZ in lKomatsu et al" ( 2008h . which gives large-scale 
temperature anisotropy IST /T — —TZ/b. An anti-correlated spectrum with (3 
large scales. 



-1 gives a positive on 



Limits are found for a_i, corresponding to anti-correlated models with (3 — —1. This could correspond 
to a curvaton scenario in which inflation is driven by a n infiaton field, but C DM perturbations are gener- 
ated by the decay of a distinct curvaton field (see e.g., iLvth fc Wand] ^^). In this case we make the 
assumption that the spectral index of the anti-uncorrelated entropy fluctuations is equal to the adiabatic 
spectral index. We do not find evidence for curvaton entropy pertubations, finding a limit from WMAP of 
Q!_i < 0.011 (95% CL), shown in Figure fT2l and in Tabled This is half the three-year limit, and places 
strong lower limits o n the possib le density of the curvaton at its decay, in this scenario, compared to the 
total energy density (jKomatsu et al., ,2008,) . If the curvaton dominated at decay, the perturbations would 
be purely adiabatic. Where has the improvement come from? The pure entropy spectrum and the cross- 
correlated spectrum both add large scale power, so a similar degeneracy is seen with the spectral index, and 
(78 1 a-s in the case where the tensor amplitude is varied. The entropy spectrum is also out of phase with the 
adiabatic spectrum, so the improved TE measurements combine with the third peak TT spectrum to tighten 
the limits. 

We also place limits on ag, correspond ing to an uncorrelated model with [3 — Q. In this case the entropy 
spectral index is set to be scale invariant. Komatsu et al. ( 20081 ) describe how this corresponds to entropy 
perturbations created by axions, which would constitute some part of the dark matter budget. The limit is 
ao < 0.16 (95% CL), ten times higher than the anti-correlated amplitude, but still preferring pure adiabatic 
fluctuations. Without the large-scale power contribution from the anti-correlated spectrum, a much larger 
amplitude is permitted, but with the same degeneracy with the spectral index and as - This has impli c ations 



for the maximum deviation from adiabaticity of axion dark matter and photons. iKomatsu et al 



(200 



provide a discussion of the theoretical implications of these limits, and those for combined data, for both 
models considered. 
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Fig. 13. — Evidence for a non-zero effective number of neutrino species, -/Veff- Left: The marginalized 
probability distribution gives iVcff > 2.3 (95% CL) from WMAP alone. The best-fit ACDM model with 
Ncs = is a poorer fit to the data than iVcff = 3, with Ax^ = 8.2. Inferred 95% limits from big bang 
nucleosynthesis (BBN) observations are highlighted. Right: Joint two-dimensional distribution for A'cff and 
the CDM density, flch^, with five-year limits in blue, compared to three-year limits in grey. The degeneracy 
valley of constant Zeq is shown dashed, indicating that the CMB is now sensitive to the effect of neutrino 
anisotropic stress, which breaks the degeneracy. 



4.2. Composition and geometry of the Universe 

4-2.1. Number of relativistic species 

Neutrinos are expected to play an important role in the dynamics of the early universe. For standard 
parameters, they contribute about 40% of the energy density of the universe during the radiation epoch 
and about 11% of the energy density of the universe at z ~ 1100 (very close to the energy density in 
baryons). Because neutrinos contribute to the expansion of the universe and stream relativistically out of 
density fluctuations, they produce a significant i mprint on the growth rate of structure and on the structure 



of the microwave background fluctuations fe.g.. lMa fc Bertschingen (|l995l )). The amplitude of these effects 



depend upon A'cff, the number of effective neutrino species. By 'effective neutrinos species', we are counting 
any particle that is relativistic at z ~ 1000 — 3000, couples very weakly to the baryon-electron-photon 
fluid, and has very weak self-interactions. Because we know neutrinos exist, we associate 'neutrinos' with 
'light relativistic particle', but note that in the strictest sense we limit only light relativistic species, as the 
cosmological constraints are sensitive to the existence of any light species produced during the big bang or 
any additional contribution to the energy density of the universe (e.g., primordial magnetic fields). 

Measurements of the width of the Z provide very tight limits on the number of neutrino species: N^, = 
2.984±0.008 (Particle Data Book), consistent with the 3 light neutrino species in the standard model. Because 
of non-thermal effects due to the partial heating of neutrinos during the annihilations, and other small 
corrections, the effective number of species is 3.0395 ( Dicus et al. 1982 : Mangano et al. 2002h . Most analyses 
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Fig. 14. — Comparison of the CMB angular power spectrum for the best-fit ACDM models with the standard 
A'off = 3.04 neutrino species (red), and with Nes = species (blue). The lower panel shows the fractional 
difference between the two spectra when iVofi is increased from to 3.04. The NcS = model has a lower 
flrnh^ in order to fit the third peak, and a lower spectral index, Us, compared to the iVoff = 3.04 model. 



of the number of neutrino species with three-year WMAP data 


Soergel et al. 


2007 


Mangano et al.ll2007; 


Hamann et al.ll2007; 


de Bernardis et al. 


2008 


) relied on combinin 



with probes of the growth rate of structure. Since one of the signatures of the number of neutrino species is 
a change in the growth rate of structure, there are degeneracies between the properties of the neutrinos and 
of t he dark energy. Neut r inos, however, leave a distinctive signature directly on the CMB power spectrum 



(see iBashinskv fc SeliakI (j2004l ) for detailed discussion): the neutrinos not only suppress the CMB peak 
heights, they also shift the acoustic peak positions. While the effects that depend on shifts in the epoch 
of matter/radiation equality are degenerate with changes in the matter density, the effects of neutrino free- 
streaming are distinct. Changes in the baryon/matter ratio and the baryon/photon ratio also have their 
own imprints on the Silk damping scale and on the acoustic scale. With five years of data, we are now able 
to see evidence of the effects of the neutrinos on the CMB power spectrum. 

Figure 16 shows the limits on the number density of neutrinos and the density in dark matter. The 
degeneracy valley, shown in the right panel, corresponds to a constant ratio of matter density to radiation 
density, or equivalently a measurement of the expansion factor at matter radiation equality: 



l + Ze 



Pc + Pb 
P7 + Pp 



40500 



1 + 0.23iVcff 



(22) 
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Parameter 


Limits 


iVeff 


> 2.3 (95% CL) 




< 1.3 eV (95% CL) 


Yp 


< 0.45 (95% CL) 



Table 6; Constraints on neutrino properties and the primordial helium fraction. 



With only 3 years of data and a lack of precision on the third peak position and height, WMAP was not able 
to make a clear detection of neutrinos (or rel ativistic species); howeve r, the data did provide a ^ 2cr hint 
of the effects of neutrino anisotropic stresses ( Melchiorri fc Serra 20061 ). Figure 16 shows that the five year 
data alone, we now constrain the number density of relativistic species: N^s > 2.3 (95% CL). By bounding 
TVoff < 10, and choosing a uniform prior on iVoff, this level of significance depends somewhat on the prior. 
We therefore test the significance of the constraint by comparing two ACDM models: one with iVoff = 0, 
and one with the standard iVcfr — 3.04. We find that the data prefer A^off = 3.04. The best-fit model has 
A(— 21nL) = 8.2 less than the Ncft — best-fit model, corresponding to evidence for relativistic species 
at > 99.5% confidence. The CMB power spectra corresponding to these two models, and their fractional 
difference, are shown in Figure fT4l The model with no neutrinos has a lower matter density, flmh^, in order 
to keep Zeq fixed. The improvement in likelihood between the two models comes from both the low-^ and 
high-^ TT spectrum, with a small contribution from the TE spectrum. We also check that this evidence does 
not go away if we relax the assumption of a power-law spectral index, by testing a model with a variable 
running, dns/dlnk ^ 0. 

Komatsu et al. ( 20081 ) combine WMAP data with other distance indicators (which constrain flch^) and 
finds a stronger limit on the number density of neutrino species: N^g = 4.4 ± 1.5. These limits will continue 
to improve as CMB measurements of the higher peaks improve. The CMB constraints on the number of 
relativistic species at redshift ~ 1000-3000 complement constraints from big bang nucleosynthesis and from 
particle accelerators. Measurements of the abundance of H elium are sensitive to the expansion rate of the 
universe during its first few minutes (jSteigman et al.lll977r ). The agree ment between the bes t fit value from 
big bang nucleosynthesis, A^cff = 3.24 ± 1.2 (95% confidence interval) (jCvburt et all |2005h . Particle Data 
Book 2007), with the best fit CMB value is another consistency check for standard cosmology. 



4-2.2. Neutrino mass 



Cosmological data places limits o n the mass of neu1: r inos. Atmospheric and solar neutrino experiments 
show that neutrinos are massive (see iMohapatra et al.l ([20051)), and measure the difference between the 
square of their masses, m^j — m^^-. Cosmological measurements constrain the sum of the masses J2 ''^f due 
to their effect o n the propagation of perturbations, on the clus tering of matter, and on the expansion rate 
of the universe teond &: Szalavlll983l: lMalll996 : IHu et aDll998h . The mass has a large effect on the matter 
power spectrum, as massive neutrinos do not cluster as well as cold dark matter, leading to a suppression 
in power on small scales. Neutrinos also affect the CMB at earlier times: if the fraction of dark matter 
that is warm is raised, acoustic oscillations in the photon-baryon plasma are less strongly damped for modes 
that entered the horizon while the neutrinos were relativistic, raising the acoustic peak amplitudes. The 
radiation-like behavior at early times also changes the expansion rate, shifting the peak positions. 



These effects are somewhat degenerate with other parameters, so CMB data alone cannot limit the 
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Xm, (eV) Zm, (eV) 

Fig. 15. — Limits on the sum of neutrino masses with the WMAP five-year data. Left: The marginahzed 
one-dimensional hmit from WMAP alone is < 1-3 eV (95% CL). This is raised by <10% with 

marginalization over a running spectral index, tensor fluctuations, or a dark energy equation of state w. 
Right: The neutrino mass is anti-correlated with as, the amplitude of matter fluctuations. 



mass as well as when combined wit h other data. With the three-year WMAP data alone the limits were 
^TTijy < 1.8 eV (jSpergel et al.ll2007l ). and < 0.66 eV when combined with other data. Since the three- year 
WMAP analysis there have been many studies of the constraints, as discussed in Komatsu et al. ( 2008h . 



The five-year WMAP data now gives an upper limit on the total mass to be ^ < 1.3 eV (95% CL), 
shown in Table [HI We have checked that this upper limit is robust to the choice of cosmological models. 
The upper limit is raised by < 10% when we include tensor fluctuations, a running spectral index, or a 
constant w ^ —I equation of state of dark energy, as sho wn in Figure ITSl Thi s dependence on additiona l 
parameters is consistent with earlier investigations by e.g.. ICrottv et al.l (|2004l ): IZunckel fc Ferreiral (|2007r ). 
A larger neutrino mass raises the amplitude of the higher acoustic peaks, hence the observed degeneracy 
with CTg (Figure [15]). Stronger constraints co me from combin i ng th e CMB data with probes of the expansion 
rate and clustering of matter at later times: Komatsu et al. ( 2008h find '^rrii, < 0.61 (95% CL) for WMAP 
combined with additional data. 



4-2.S. Primordial Helium Abundance 

In most cosmological analyses the primordial helium abundance is fixed to be Yp — 0.24, motivated by 
observations discussed in Sec 13.2.91 The effect of the abundance on the CMB spectrum is small, but provides 
an independent cross-check of the BEN results, and probes for any difference between the helium abundance 
during the first few minutes, and after 300,000 years. The abundance affects the CMB at small scales due 
to the recombination process. The number density of electrons before recombination depends on the helium 
fraction through rig = nf,(l — Yp) where ni, is the baryon number density. Changing the electron number 
density changes the mean free path of Compton scattering, which affects the Silk damping scale. A larger 



- 35 - 



Yp increasingly damps the power on small scales, as shown in iTrotta fc Hansen! (j2004r ). 



Co nstraints from the firs t -year WMAP data were presented in iTrotta fc HansenI (120041): iHuey et a,l 



(|2004l ): llchikawa fc Takahashil (|2006l ). with 99% upper hmits of Yp < 0.65 inferred (|Trotta fc Hansenll2004f) . 
A subsequent analysis of the three- year data gave Yp < 0.61 at 95% confidence, tightened to 0.25 ±0.10 with 
small-scale CMB data (jlchikawa et ahlbooj . We now find Yp < 0.45 (95% CL) with the five-year WMAP 
data. Higher values allowed by the thre-year data are disfavored with a better measure of the third acoustic 
peak height. With future small-scal e CMB measurements, for example fr om the Planck satellite, constraints 
should significantly improve (|Ichikawa et al.ll2007t iHamann et al.ll2008b( ). 



4-2. 4- Curvature of the universe 



In combi nation with other d ata, WMAP observations place strong constraints on the geometry of 
the universe (jSpergel et al.l 120071 ). The CMB measures with high accuracy the angular scale at which 
acoustic oscillations are imprinted at the last scattering surface, 6** = 0.5952 ± 0.0017 °. However, this 
alone does not provide a good measure of the geometry, as there is a degeneracy with the expansion rate 
of the universe since last scattering. This is shown in Figure I16[ indicating the degeneracy between the 
dark energy density ^Ia and the curvature ilk- With WMAP alone the curvature is weakly constrained, with 
marginahzed limits = -0.099l[!;?gg, and Ha < 0.76 (95% CL), assuming a Hubble prior of 20 < Hq < 100 
and CI A > 0. Without this prior on the positivity of CIa, limits on the curvature are weakened. The 
same degeneracy is seen, although slightly broadened, when the dark energy equation of state is allowed to 
vary. However, in both cases the Hubble constant decreases with increasingly negative curvature, taking 
values inconsi stent with observatio n. This degeneracy can be used to con strain the curvature by combining 
observations ( Jungman et al. 1996h . In the three- year WMAP analysis, Spergel et al. (2007) showed that 
the degeneracy is truncated with the addition of only one piece of additional cosmological data (Type la 
supernovae, or the HST meas urement of the Hubble constant, or galaxy power spectra), tightly constraining 
any deviations from flatness. Komatsu et al. ( 20081 ) draw similar conclusions with currently available data, 
and discuss the current limits on the spatial curvature from recent observations. 



4-2.5. Dark energy properties 



The ACDM model requires a non-zero dark energy density Ua = 0.742 ± 0.030 to fit the data, which is 
assumed to be in the form of a cosmological constant. We do not have a n explanation for th is component 



of the universe. A natural explanation could be a vacuum energy density ( Carroll et al. 1992h . but if so, we 



are faced with the fine-tuning problem to explain its observed value, 120 orders of r nagnitude smaller than 
expected from field theory arguments. A lternative explanations incl ude quintessence ([Peebles fc Ratr al ll988l: 



Wetterich 1988; Ferreira fc Jovcel 19981 ) or modifications to gravity ( Deffavet et al. 2002 ). Testing the dark 
energy equation of state today, and as a function of cosmic time will help identify the possible explanation. 

The CMB by itself cannot place strong limits on the equation of state w — p/ p, but by measuring the 
acoustic peak positions and heights, and constraining flmh^ with the third peak, limits the range of models 
to a degeneracy between ilm and w, shown in Figure 1161 The dark energy in these models is allowed to 
cluster. With a prior on the Hubble constant Hq < 100, WMAP alone places weak limits w — —1.06 



0.41 
0.42J 



with Qa = 0.73_q']^]^. If flatness is not assumed, the WMAP data cannot constrain w or Qa due to the 
geometric degeneracy, also shown in Figure fTBl However, the situation is significantly improved when WMAP 
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Fig. 16. — Left: The points show the set of non-flat models consistent with the WMAP data, colored by the 
Hubble constant values. WMAP measures the acoustic peak scale to high accuracy, but does not constrain 
the curvature, fifc, by itself. However, the highly curved models have a low Hubble constant, inconsistent 
with observation. Right: Constraints on the dark energy equation of state, w, and the dark energy density, 
Ha, from WMAP alone. With a Hubble constant Hq < 100, weak limits can be placed on it; in a flat universe, 
shown by the blue contours, but the dark energy density and equation of state are unconstrained (with the 
95% confidence level shaded grey) if the assump tion of flatness is rela xed. Limits are signiflcantly improved 
when WMAP is combined with additional data ([Komatsu et al.ll2008| ). 



is combined wi t h astr onomical data measuring the expansion rate and clustering of matter at late times. 
Komatsu et al. ( 2008[) discuss limits obtained from various data in combination with WMAP, and find w 
constrained to be —1 to within 6% for a flat universe and constant equation of state. 



5. Conclusions 

The simple six parameter ACDM model continues to fit the WMAP data. With five years of observations, 
we have better measured both the temperature and polarization anisotropy of the CMB. This has allowed 
us to measure with smaller errors, compared to the three-year analysis, the third acoustic peak in the 
temperature spectrum, and the low-£ polarization signal, leading to improved constraints on the cosniological 
parameters describing the contents of the universe, and the primordial fluctuations that seeded structure. 
The observations continue to be well fit by the predictions of the simplest inflationary models, with a scale- 
invariant spectrum of fluctuations disfavored. Consistency with the TE cross-correlation spectrum, now 
measured with better accuracy, provides additional confidence in this simple model. 

We have detected the optical depth to reionization with high significance. This measurement implies 
that reionization of the universe likely took place gradually, as it constrains a sudden reionization to be earlier 
than consistent with other observations. With more data, it will become possible to use the polarization data 
to better quantify the ionization history. Given the improvement in this measurement, and with a view to 
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interpreting future large-scale polarization measurements, we develop an alternative way to remove Galactic 
foregrounds from low resolution polarization maps, which includes marginalization over uncertainties in the 
Galactic signal. We find consistent results using this method and the standard template-cleaning method. 

Considering a range of extended models, we continue to find that the standard ACDM model is con- 
sistently preferred by the data. The improved measurement of the third peak now requires the existence 
of light relativistic species, assumed to be neutrinos, at high confidence. The standard scenario has three 
neutrino species, but the three-year WMAP data could not rule out models with none. The ACDM model 
also continues to succeed in fitting a substantial array of other observations. Certain tensions between other 
observations and those of WMAP, such as the amplitude of matter fluctuations measured by weak lensing 
surveys and using the Ly-a forest, and the primordial lithium abundance, have cither been resolved with 
improved understanding of systematics, or show some promise of being explained by recent observations. 
With further WMAP observations we will better probe both the universe at a range of epochs, measuring 
fluctuation characteristics to probe the initial inflationary process, or other non-inflationary scenario, im- 
proving measurements of the composition of the universe at the recombination era, and characterizing the 
reionization process in the universe. 
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We have switched from the pixel likelihood to a Blackwell-Rao (BR) estimator at low ell, so we perform 
several tests to verify that the new likelihood 1) is consistent with the pixel likelihood, 2) is insensitive to 
several choices of input data, and 3) is properly converged. Many of our checks involve comparing likelihoods. 
Because the likelihood is a function on a high dimensional space, we do not check it everywhere. Instead 
we choose some fiducial spectrum and compute slices through the likelihood. We compare these conditional 
likelihoods, typically in the form of maximum likelihood values and full-width-half-max error bars between 



A. 



Low-^ TT likelihood cross checks 



resolution parameter 
smoothing FWHM 
concise /pixel [jjK] 
^max sampled 
^max conditioned 



4 

9.1831° 



1.0 

32 

48 



5 
5.0° 
2.0 
51 
96 



Table 7: Parameters used for sampling the \ow-i TT likelihood at two different resolutions. 
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Fig. 17. — The angular power spectra of signal and noise components in the smoothed ILC map used for 
Gibbs sampling the low-^ temperature spectrum. Uncorrelated noise, at 2 iiK per pixel (red), is added to 
the smoothed ILC map to speed up the sampling, and then is assumed to be the only noise present. This 
assumption is inaccurate at low as it ignores the true noise (blue), but the error is negligible since it is 
significantly lower than cosmic variance (green). 

different likelihood methods. We also compare ACDM cosmological parameters obtained using the 

different inputs. 

The first comparison is between the pixel-based code and the Gibbs code. Figure [1] showed a compar- 
ison of the Ci likelihoods, with no significant differences, even though they use different resolutions. The 
cosmological parameters are unchanged for the two codes. The next comparison is between likelihoods from 
Gibbs sampling at two different resolutions, iVgide = 16 and 32. Table [7] lists details of runs at these two 
different resolutions, including the standard deviation of Gaussian white noise added to the smoothed maps. 
The power spectrum of the white noise for the A^sjdc = 32 case, is compared to the input noise and cosmic 
variance in Figure [T71 Power spectra are sampled up to ^max (sampled), and the sampled CMB skies are 
constructed with the sampled power spectra and by conditioning on a constant fiducial spectrum out to 
a higher value of €,„ax (conditioned) . In Figure [18] the estimated spectra are compared, using the fiducial 
KQ85 mask. There are some small differences, but the spectra are consistent. 

We also show in Figure [TS] that the spectra obtained are almost identical for two different realizations for 
the added uncorrelated white noise, and that using the larger KQ80 mask, compared to the fiducial KQ85, 
has a small effect on the spectra and cosmological parameters, consistent with noise. Finally we compare 
results for different ^max used in the BR estimator. The standard set-up uses fmax = 32; using the low-£ 
likelihood up to £ = 51 has almost no effect on cosmological parameters. 
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Fig. 18. — The level of variation in the low-ell TT Gibbs likelihood that arises from different input parameters. 
The spectrum estimated from the standard input (black) is compared to results obtained using a different 
white noise realization (blue), using a larger mask (KQ80,red), and degrading to A^side = 16 (green) rather 
then 32. The likelihood does fluctuate with these changes, but has a negligible effect on cosmological 
parameters. The values shown at each multipole correspond to maximum likelihood values obtained by 
fixing the spectrum of other multipoles at fiducial values. The error bars show where the likelihood is at 
50% and 5% of its peak value. 



B. Parameter estimation details 



B.l. Sampling method 



In the primary pipeline, only one chain is run for each model, rather than the typical four or eight 
parallel chains used in most MCMC cosmological analyses. This is possible as our spectral convergence test 
only requires a single chain. The starting points are picked as good-fitting points from previous analyses 
(e.g., a good-fitting value from the three-year WMAP analysis), or previous test chains, which means that no 
initial burn-in need be removed. Starting afresh with the five-year data, one can always use a point lying in 
the WMAP three-year limits as a starting point. This would not be the case for entirely new distributions, 
or much improved data, in which case we run a short initial exploratory chain to find the high likelihood 
region. 

The covariance matrix from which trial steps are drawn, is chosen using a best guess for the covariance 
C of the distribution being sampled. For a multi-variate Gaussian with covariance C, the optimal trial 
covariance matrix, if all parameters are sampled simultaneously, is = (ctt/o'o)^CI, where 



(Bl) 



CTT _ 2.4 

(Jo ^ Vd 

for D-dimcnsional distributions ( Gelman et al. 1996 : Hanson fc Cunningham 1998t Dunklev et al. 2005h . 
This relation also holds for somewhat non-Gaussian distributions. At each step we draw a vector G of D 
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Gaussian unit variance zero mean random variates, and compute a trial set of parameters x^, starting from 
the current position x^, where 

XT = Xi + VCtG. (B2) 

Using an appropriate covariance matrix can speed up the samphng time by a factor of hundreds. In practice 
we have a good idea of the covariance from previous cosmological analyses, so this is always used as a starting 
point. With significantly different data, or with new parameters, a best guess is made. The chains are then 
run for a few thousand steps, and then updated if the matrix is inadequate. This is determined by the 
acceptance rate of the chain, which should be ^ 15 — 25% and by the chain efficiency, using the spectral 
test described in the next section. A second or third update may be required for models such as those with 
curved geometries with variable dark energy equation of state. Once a good covariance matrix is found, the 
chains are run for typically 20,000 steps, and then tested for convergence. Convergence for ACDM chains 
typically takes place after only ~6000 iterations, although at this point the ID and 2D distributions are 
noisy. In practice we then run all the chains for longer than the convergence limit, typically for 100,000 
iterations, in order to get well sampled histograms. The chains are not thinned before analysis. 



B.2. Spectral convergence test 



A chain has 'converged' when its statistical properties reflect those of the underlying distribution with 
sufficient accu r acy. To determine this stopping point, we use the spectral convergence test described in 
Dunklev et al.l (j2005l ) in our main parameter pipeline. The power spectrum of each parameter of the chain 
is used as a diagnostic, to check whether the chain has (1) sampled the distribution in such a way that it 
is unbiased by correlations, and (2) s ampled enough points that statistics can be estimated with sufficient 
accuracy. The Gelman & Rubin test ( Gelman fc RubinI 1992 ). commonly used in cosmology, can sometimes 
fail to test the first point, producing a false positive. 



To estimate the power spectrum P{k) from a chain of length we construct Pj = |ajp 
27rfc/7V, by taking the discrete Fourier Transform of the chain of values for each parameter, x, 

x{n) exp[i2TT{jn/N) 

71=0 



where j 



(B3) 



where — (A^/2 — 1) < j < N/2. x{n) is the value at each iteration n, so chains stored in weighted format 
are converted to unweighted arrays for analysis. Since the Metropolis algorithm produces chains which are 
correlated on small scales, the power spectrum tends to a white noise spectrum on large scales, and turns over 



to a spect rum with suppr e ssed p ower at large fc, with the turnover position reflecting the inverse correlation 
length. In lDunklev et al.l ( 20051 ) it is shown that the spectrum can be fit by the following template: 



P(fc) = Po 



(fcVfc)'^ 



(fc*/fc)" + 1' 



(B4) 



with Pq giving the amplitude of the white noise spectrum in the fc ^ limit, fc* indicates the position 
of the tur nover to a different po wer law behavior, characterized by a, at large fc. This model is shown in 
Figure 3 of bunklev eTZI |2005[ l. fitting the noisy spectrum of a parameter from a chain. Th e model fits the 



noise- averaged spectrum of a real chain obtained from Monte-Carlo simulations, also shown in lDunklev et al 
(l2005l) . 
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To fit the parameters In Pq , k* and a to P, using least squares for a finite chain, we have 



InP, = ln[Po] +ln 



(7Vfc*/27rj)" 



1 + {Nk* /27r jy 



7 + r 



1 ' 



(B5) 



where 7 = 0.577216 is the Euler-Mascheroni constant, and rj are random measurement errors. The parame- 
ters are fit over the range of Fourier modes 1 < j < 10 j*, for a spectrum with j* — k*{N/2TT), so we iterate 
twice to converge on the j* limit. 

For convergence, the largest scales probed must be in the white noise regime P{k) fc", defined by 
the requirement j* > 20 for each parameter. This insures that the correlated points are not biasing the 
distribution and indicates that the chain is drawing points throughout the full region of high probability. To 
test sufficient accuracy, we require the convergence ratio r — a^/aQ to be less than 1% for each parameter, 
where a'^ is the variance in the sample mean, and (Tq is the variance of the sampled parameter. The Gelman 
& Rubin test incorporates a similar ratio: their R statistic r oughly translat e s to R ~ 1 + r, but the quantity 
is calculated using multiple parallel chains. It is shown in IPunklev et al.l ( 2005 ) that r can be estimated 
using a single chain, since estimating the sample mean variance of a long chain, with zero mean, is equivalent 
to estimating P(fc) at fc = 0: 



1 

N 



P{k = 0). 



(B6) 



In practice we rescale each parameter to have zero mean and unit variance before computing its power 
spectrum. Then we estimate the value of Pq for each parameter, to compute r = Pq/N. We require r < 0.01, 
but in practice obtain much smaller values, typically with r < 0.001. 
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